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ABSTRACT 


The pyrolysis of pure and fire retarded bulk cellulose samples in a 
nitrogen atmosphere is studied. Sodium hydroxide is used as the fire 
retardant. The study is directed toward determining the effects of the solid 
phase fire retardant on the burning of cellulose. The material property of 
pure and fire retarded cellulose which most directly affects its burning 
behavior, the heat of gasification, is measured by using a specially designed 
pyrolyzing chamber. Theoretical results are obtained for a one-dimensional 
pyrolysis wave propagating into cellulose by using a  finite-difference 
calculation. The experimental results are in qualitative agreement with the 
computations, but give a lower value for the heat of gasification. The mass 
fraction of noncombustible volatiles in the total volatile products of 
pyrolysis, the stoichiometric oxygen/fucl ratio of the combustible volatiles 
and the heat of combustion per unit mass of combustible volatile are also 
obtained as experimental results. These data have been obtained because they 
are needed as inputs in diffusion flame calculations. The experimental data 
show that sodium hydroxide acts as a catalyst in the pyrolysis of cellulose. 
Its addition leads to decreases in the heat of gasification and the mass 
fraction of noncombustible volatiles in the total volatiles, and to increases 
in the stoichiometric ratio and the heat of combustion of combustible 
volatiles, thereby has a dual effect on cellulose burning. 

The experimental techniques developed for the study will be useful in 


measuring the combustion properties of other charring materials and will 


-Vil- 


assist in the development of scientifically based test methods. The 
experimental results are useful in determining the physical mechanisms by 
which solid phase fire retardants act to suppress flaming combustion, and 
will help to provide a rational basis for the design of materials with 


improved fire safety characteristics. 
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ie Heat of reaction 
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CHAPTER 1 


INTRODUCTION 


1.1 General Background on Cellulose Pyrolysis 


Cellulose is a major constituent of combustible materials available as 
fuel in fires, for example, Fremy (1868) showed that wood is 40% cellulose. 
The combustion of cellulose, therefore, has always been closely associated 
with the development of mankind. An essential feature of the combustion of 
cellulose is its pyrolysis, or thermal degradation. Considerable effort has 
been expended in the study of cellulose pyrolysis because of its central 
importance in fire research and cellulosic materials utilization. The study 
of cellulose pyrolysis is closely associated with practical problems such as 
textile flame-retardancy, energy and chemical production from. thermal 
conversion of biomass and air pollution research. 

A cellulose molecule consists of long chains of repeating glucosan 
units, each having a formula of C6H1905.- Cellulose seems to be divided 
up into periodically spaced crystalline or highly ordered regions, in which 
the molecules are linked together by hydrogen bonds, and into amorphous or 
highly disordered regions, where this kind of bonding is mostly absent. 

The literature on cellulose is mainly concerned with its thermal 
degradation, although some burning studies have also been reported. An 
extensive survey of this literature is given by Shafizadeh (1968); more 
recent work is discussed in Lewin and Basch (1978). 

Pyrolysis encompasses the processes by which gaseous fuel is liberated 


to support fire by the breakdown of the fuel constituents under the 


ae 


influence of heat. The chemistry of cellulose pyrolysis is very complex. 
Different mechanisms have been suggested by numerous investigators. The 
One pointed out by Kilzer and Broido (1965) has been generally accepted. It 
maintains that pyrolysis of cellulose proceeds by two distinct competing 
paths, which involving at least three concurrent and consecutive processes, 
each including one or more reactions. As illustrated, in Fig..lsh) the two 
paths include: (1) the dehydration and char-forming reactions, which are 
initiated at lower temperatures and continue at elevated temperature; and (2) 
depolymerization and levoglucosan or tar-forming reactions, which gather 
momentum at temperatures above 250°C and result in rapid volatilization. 
The first process, which is slightly endothermic, is attributed to an 
inter-molecular reaction with loss of water to form a "“dehydrocellulose”. 
During the second endothermic process, which competes with the first 
process for native cellulose, levoglucosan, the major constituent of the tar, is 
postulated to arise from the rearrangement of 1, 4-anhydro-a-D-glucopyranose 
formed in an "unzipping" reaction. The tar then breaks down into highly 
flammable volatiles, which support flaming combustion, and into. char. 
Finally, emthes: ¢xothernmicusprocess,»:1S,.,..ascribed,. tom,reactions,. of patie 
"dehydrocellulose", which experiences carbon-carbon and carbon-oxygen bond 
ruptures and hydride ion transfers to produce volatile carbon-containing 
compounds and hydrogen, and inter-molecular condensations to produce 
char. 

Although the nature and extent of many individual reactions involved 
in cellulose pyrolysis are not known, or sufficiently defined, it has been 


found by Lipska et al. (1966), McCarter (1972), Lee et al. (1977 and 1981), 
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Kashiwagi et al. (1981), Cullis et al. (1983), Mok et al. (1983), Nunn et 
al. (1985) and Ohlemiller (1987) that these reactions are highly influenced 
by temperature, external radiant flux, external heating rate, total time of 
heating, sample type, amibient oxygen concentration, water content, and 
pressure. Lipska and Parker (1966) studied the isothermal pyrolysis of 
cellulose in nitrogen over a temperature range of 250 - S00 G: and found 
that decomposition and volatilization occuring during pyrolysis, increases 
with increasing temperature. McCarter (1972) performed pyrolysis studies of 
cellulose at rates approaching those in burning, and found that the 
activation energy at high temperature differs from that previously reported 
at lower temperature. Mok et al. (1983) investigated the effects of pressure 
on cellulose pyrolysis and found that increased pressure favors char 
formation. 

Fire retardation is essentially an interruption or suppression of the 
burning process. Research on the fire retardation of cellulose has had a 
long history with the aim of preventing destructive fires. Tneresise a 
voluminous literature on the effects of fire-retarding additives on _ the 
thermal degradation of polymeric materials. An extensive discussion of fire 
retardancy is given in the Encyclopedia of Polymer Science and Technology 
(1967). Systematic coverage of fire retardants is presented in Fire Safety 
Aspects of Polymeric Materials, Vol. 1. (1977). Current commercial 
fire-retarded materials have been reviewed by Bauer (1978). The 
previously mentioned reviews by Shafizadeh (1968) and Lewin and Basch 
(1978) also discuss effects of fire retardents on cellulosic materials. 


The theories, which have been advanced to explain the effect of fire 


vis 


retardants on cellulose are categorized by Shafizadeh (1968) and Mark et al. 
(1975) as follows. 

(1) Coating therories. These attribute the retardant action to a coating 
of the fibers by. the melted or foamed retardant, or both; restrictingsen. 
escape of volatile products of cellulose pyrolysis and the access =sor 
atmospheric oxygen to the reaction zone. 

(2) Thermal theories. These attribute the retardant action to a lowering 
of cellulose temperature, brought about by: (a) insulating cellulose by a 
coating of melted or foamed chemical, or both; or (b) absorbing the heat, 
which would otherwise reach cellulose, by endothermic reaction in the 
retardant; or (c) increasing thermal conductivity of cellulose and thus 
slowing down the temperature rise. 

(3) Chemical theories. These attribute the retardant actions a1eee 
modification of the initial pyrolytic reactions in cellulose pyrolysis through: 
(a) catalysis of decomposition reactions leading to dehydration and char 
formation, thus increasing the production of CO,5, HO and char; and (b) 
inhibition of decomposition reactions leading to tar formation, and thus 
decreasing the production of flammable volatiles. 

(4) Gas phase theories. These attribute the retardant action to 
modification of the volatiles from cellulose pyrolysis by gases evolved 
during decomposition of the retardant, such as (a) diluting the flammable 
volatile decomposition products with nonflammable gases, such as CO, and 
HO; and (b) introducing decomposition products from the retardant that act 
as chain breakers to terminate free-radical chain reactions in the flame. 


Many researchers, such as Kilzer et al. (1965), have believed that the 
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two competitive paths feature of pyrolysis, as shown in Fig. 1.1., applies to 
fire-retarded cellulose as well, but McCarter (1972) believes that the two 
main paths are sequential, not competitive, for retarded cellulose. Tang et 
al. (1964), Kilzer et al. (1965), McCarter (1972) and Shafizadeh (1968 and 
1975) believe that fire retardants work by catalysis of the dehydration path, 
thus increasing the production of char, water and CO5 at the expense of 
tar, which releases the flammable volatiles. 

A notable early study on the effect of flame retardants on cellulose 
was performed by Tang and Neill (1964). They found an increase in the 
yield of char, and a decrease in the heat of pyrolysis when fire retardants 
are added. In another study of cellulose pyrolysis, Lipska and Martin (1971) 
found strong effects of acid treatment on the isothermal pyrolysis rate. 
These types of studies are continuing; a recent example is Rogers and 
Ohlemiller (1981) who investigated the degradation of cellulosic material of 
the type used in home insulation. Paker and Lipska (1969) and Mok et al. 
(1983) believe that some of the char is produced by the action of retardants 
on levoglucosan. Drews and Barker (1974) found that the chemical nature of 
the char produced by pyrolyzing cellulose was independent of the different 
types of organophosphorous compounds added as retardants. Woodley (1971) 
found that a retardant such as KHCO3 causes changes in pyrolysis vapor 
products of cellulose. The pyrolysis products depend on the concentration, 
rather than the type of retardants. Kilzer et al. (1965) speculated that 
water helps cool down the condensed phase to a temperature at which it no 
longer supplies flammable volatiles. Mok et al. (1983) believe that 


retardants tend, overall, to lower the decomposition temperature of cellulose, 
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and thereby modify its pyrolysis so that char residues are increased, and the 
amount and/or flammability of the evolved vapors are decreased. 

A number of fire retardants, primarily inorganic salts, have been 
studied in the literature cited Rate. The most comprehensive comparative 
study of fire retardant effects on cellulose was made by Shafizadeh et al. 
(1975) -who*ranked) 2], sretardants. on the » basis».of, theirs,effiect ,ong seus 
production and on the heat of combustion of the gaseous pyrolysate. The 
retardants which span the range of effectiveness given by Shafizadeh, in 
order of increasing effectiveness, are (NH4)> CO3, NaHCO3, Na,CO3 and 
NaOH. 

Lipska and Martin (1971) studied retarded cellulose in a _ nitrogen 
atmosphere> and;),..found , that, the. presence,, of, oxygen, ,even Anaetow 
concentrations, will increase the rate of decompostion in the initial period 
of pyrolysis. They also found that the overall pattern of decomposition of 
pyrolyzed cellulose treated with basic and neutral additives is similar to 
that of the untreated cellulose, but the pyrolytic decomposition of cellulose 
treated with acidic retardants takes on different characteristics. They 
suggested’ “thatver the §hrextremely |, rapid «’.decomposition...prates., «of awsthe 
acidic-salt-treated samples may be caused by the combined efforts of 
hydrolysis and pyrolysis. 

As a result of rapidly developing computational. techniques, theoretical 
modelling of cellulose pyrolysis has become popular. Kanury and Blackshear 
(1970) proposed a simple linear model to compute the thermal properties of 
partially charred, woody solids, given information A thermal conductivity 


and specific heat as a function of the extent of charring. Roberts (1971) 
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conducted a thermal analysis of wood pyrolysis and concluded that the 
conventional assumption of constant thermal diffusivity and no _ internal 
convective heat transfer cannot be justified for the entire wood pyrolysis 
processes. In the model developed by Kung (1972) for wood pyrolysis, 
thermal decomposition was assumed to occur by a first order reaction with 
an Arrhenius temperature dependence. A parametric study was made of the 
sensitivity of such results as temperature-dependent mass loss rate to changes 
in the assumed values of the kinetic constants and heat of gasification. This 
type of amalysis was extended by Kansa et al. (1977) to include a 
momentum equation for the motion of the pyrolysis gases in the interior of 
the solid, thus allowing the build-up of internal pressure to be calculated. 
Even though numerous models have been developed to study cellulose 
pyrolysis and predict its results, at present there is no model which 
accurately describes the entire cellulose pyrolysis process. However, most 
researchers agree that the thermal properties of cellulose and char, and the 
"heats" of pyrolysis are the factors which most strongly affect cellulose 


pyrolysis, and thus need to be studied further. 


12 + Previous Work on “Heats” of Pyrolysis 


In the literature, several kinds of "Heat" terms have been used to 
describe the energy changes occuring during cellulose pyrolysis. These are: 
(1) the heat of pyrolysis, hp which is defined as the difference in enthalpy 
between the pyrolysis products (volatiles plus char) at the cellulose sample 


surface temperature and the original cellulose at ambient temperature; (2) 


Pg i 


the heat of reaction, h,, which is defined as the difference in enthalpy 
between the products and reactants at ambient temperature; and (3) the heat 
of gasification, hg, which is defined as the ratio of the net surface heat 


flux to the volatile generation rate due to pyrolysis. Whereas Ay and h, 


are thermodynamics quantities, hg ismnot While hy and h, can be 
calculated from other thermodynamic quantitites, hy must be obtained by 
measurement. By conventions, Ap: h, and hg aré’ /ipositive’ iOraean 


endothermic reaction. 

Experimentally, DTA (Differential Thermal Analysis) and DSC 
(Differential Scanning Colorimetry) are usually used to obtain the heat of 
pyrolysis and the heat of reaction. Tang and Neill (1964) used DTA to 
measure the heat of reaction for both unretarded and retarded cellulose 
samples. They found that retarded samples have lower heats of reaction than 
unretarded samples, and give a value of 368415 J/g (endothermic). They 
speculated that this is the result of the endothermic depolymerization and 
dehydration of cellulose and the exothermic polymerization and 
aromatization of pyrolysis products. Roberts (1971) conducted a survey on 
the heat of reaction during the pyrolysis of cellulose. The data reported, 
which are from DTA and other experiments, cover a range from 370 J/g 
(endothermic) to -80 J/g (exothermic). It was concluded that the width of 
the range of results is due to variations in experimental method and in the 
conditions of pyrolysis. Mok et al. (1983) used DSC to investigate the effects 
of pressure and purge gas flow rate on the heat of cellulose pyrolysis. It 
was found that the heat of pyrolysis decreases with increasing pressure and 


decreasing flow rate, and that preheating of cellulose samples produces a 
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lower heat of pyrolysis. 

Lee et al. (1977) performs experiments on wood pyrolysis at fire-level 
surface heat fluxes using laser simulation, and found that the overall heat 
of reaction varies from endothermic in the case of low heating rate to 
exothermic for highmiheatmemsrate, ranging trom® 610 J/g to’ -1714 "J/g, 
Tewarson (1980) obtained experimental data on the heat of gasification for 
oak which gave 1700 J/g for external heat flux ¢ 46 kW/m2, and 5500 J/g 
for external heat flux > 46 kW/m?. 

Kung et al. (1973) developed a mathematical model for the pyrolysis 
of a wooden cylinder, which, in addition to the usual consideration of 
transient heat conduction and Arrhenius reaction rate, included terms 
allowing for internal heat convection of volatiles and variable thermal 
properties. The model was used together with experimental values from 
Roberts and Clough’s (1963) work to yield a heat of reaction of 202 J/g 
(endothermic). Calculated results from the model also show that the pyrolysis 
rate increases with an assumed increase of heat of reaction. The heat of 
reaction is deemed as one of the important parameters in wood pyrolysis. 
Sibulkin (1986) examined the relationships between the heats of reaction, 
pyrolysis and gasification for vaporizing and charring materials. It was 


shown that the heat of pyrolysis, h,, is equal to the heat of gasification, 


iD? 
hg, only for steady burning of thermally thick, vaporizing material, such as 
PMMA (polymethyl methacrylate); for charring materials, such as cellulose, 


there is no simple relationship between h, and hy. Solutions for pyrolysis of 


8 
cellulose were also obtained by him using a one-dimensional numerical 


model with material properties estimated from the literature. It was found 
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that after an initial transient, he is about twice Ap: For a semi-infinite slab 
the value of he remains nearly constant; while for a finite slab the value 
of hg decreases with time to values less than hp: Changes in the assumed 
chemical reaction rate and the incident heat flux are found to haveua 
minor effect on hg. The heat of gasification is also relatively insensitive 
to changes in the char density and thermal conductivity. The variation 


of hg with time is nearly independent of the assumed heat of reaction, h,, 


but the magnitude of hg varies almost directly with the value of h,. 


1.3 Previous Work on Products of Pyrolysis 


In order to understand pyrolysis mechanisms and develop more 
effective fire retardants, many studies have been done on the pyrolysis 
products of wood and pure and retarded cellulose. Generally, pyrolysis 
products, of acelluwiosesminelude char star. CO, CO>, H,0 and other volatiles 
according to the mechanism described in Section 1.1. 

Instruments  ~suchsaas, .GC. (Gas. .Chromatograph) 2 and, VMS so (Niase 
Spectrometer) are extensively used in the experiments for pyrolysis products 
measurements. Lipska and Martin (1971) performed GC analyses of the 
pyrolysis of cellulose in nitrogen and found a complex range of products 
produced. Wodley (1971) used GC and MS to measure the pyrolysis products 
of both unretarded and retarded cellulose and levoglucosan. The results 
indicated that fire retardants decrease the number of compounds formed 
during pyrolysis. The comparison of the pyrolysis products gencrated under 


different atmospheres (Ny and air) in the temperature range of 330-400°C 


alge 


showed that the formation of pyrolysis products is independent of either 
atmosphere or temperature. The pyrolysis products formed for different fire 
retardants and concentrations indicated that the products depend more on 
the concentration than on the type of retardant. Drews and Barker (1974) 
investigated the chars produced from the burning of a series of cotton 
fabrics treated with different types of phosphorus compounds, and found 
that the chemical nature of the char formed seems to be independent of the 
chemical structure of the phosphorus compound. Kashiwagi et al. (1981) 
measured the pyrolysis products of thermally thick white pine using GC. 
The experiments were done for different radiant fluxes in both nitrogen 
and nitrogen-oxygen atmospheres. The results indicated that the yields of 
char, CO, CO, and THC (total hydocarbon) increase with increasing radiant 
flux inysawenitrogens-atmosphere. In oxygen-containing atmospheres, the 
effects of the radiant flux on the evolved gaseous products are different. 
CO, production increases significantly with increasing radiant flux, whereas 
CO evolution does not. It was pointed out that ambient oxygen also 
significantly increases char depth and changes the volatile products 
distribution during pyrolysis. It was concluded, therefore, that the effects 
of ambient oxygen on the gasification process cannot be ignored. are 
different. The CO, production increases significantly with increasing 
radiant flux, but CO evolution does not change significantly at different 
radiant fluxes. It was pointed out that ambient oxygen also significantly 
increases char depth and changes the volatile products distribution during 
pyrolysis. It was concluded, therefore, that the effects of ambient oxygen on 


the gasification process cannot be ignored. 
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Chan et al. (1985) performed experiments on product formation in the 
rapid pyrolysis of thermally thick wood particles. The results show that 
the thermal properties, heat flux, and particle size, interact in a nonlinear 
way to alter the ultimate yields, which include char, tar, water, and gases 
such as CH4, C5H4, *€O,"CO>5, and’ H50. It was also found that as 
particle size and heat flux are varied in the experiments, the gas products 
contain changing proportions of chemicals such as CHy and CH, as 
functions of time, which indicates a potential for optimization of product 
distribution from thermally thick biomass particles in the event that costly 
size reduction must be avoided. Roy et al. (1985) studied pyrolysis of 
aspen poplar under vacuum. The pyrolysis products obtained are char, water, 
gases such as CO, CO>, CHy, and pyrolytic oils which contain liquid phase 
Organic compounds, such as acetic and formic acids. Parameters such as 
temperature, particle size, heating rate, total pressure, and residence time of 
gases and vapours in the reactor were examined. Low pressure, short 
residence time and high heating rate were found to favour the pyrolytic oil 
yield, and the yield of char was correspondingly low. Nunn et al. (1985) 
studied the effect of temperature on the yields, composition and rates of 
evolution of major products from the rapid pyrolysis of sweet gum 
hardwood. It was found that the pyrolysis products include char, tar, CO, 
COs, H50,"CH4; C457; HCHO, and acetic acid. The experimental results 
indicated that tar is the major pyrolysis product at temperatures above 
800K; secondary cracking of this tar contributes significantly to the yields 
of’ CO, "CHa, Cag atid” most of the other light gaseous products at 


temperatures above 900 - 950K. On a weight basis, CO is the dominant 


so}. 


gaseous product above 850K, while CO, and water dominate the gaseous 
products below 850K. Total gas production (including water) increases 
monotonically with increasing temperature, where light hydrocarbon gases 
such as CHy, C5Hy and C3H¢ are relatively minor products. Ohlemiulervct 
al. (1987) did experiments on wood pyrolysis at fire-level heat fluxes in 
atmospheres of varying oxygen concentration. The results show that the 
major pyrolysis products, which include CO, CO, water, tar and THC vary 
rapidly ‘with time Sonlyeertnemittstiecs Or 3,..minutes but their rate of 
emission continues to vary more slowly with prolonged heat exposure. It 
was also indicated by GC and MS measurements that the organic condensate 
Or tar is extreamely complex in composition, probably exceeding 700 
constituents. 

On the other hand, theoretical work on predicting pyrolysis products is 
in steady progress as well. Chan et al. (1985) developed a one-dimensional 
model to predict volatiles release rate and composition as a function of 
biomass physical properties. The predicted shape and amplitude compare 
favorably with the experimental data. However, further refinement of 
the current kinetic mechanism and thermal property variation with 
temperature in the model is needed. Nunn et al. (1985) did a linear 
simulation of whole wood pyrolysis behavior based on kinetic data for the 
pyrolysis of cellulose and lignin. It was found that the model predicts total 
weight loss (i.e. total weight of products), quite satisfactorily, but breaks 
down for individual products and total gas yield. inerchores further model 


improvement is needed. 
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1.4 Application to Diffusion Flames 


Flames are broadly classfied as pre-mixed and diffusion flames. A 
pre-mixed flame occurs if fuel and oxidizer are homogeneouly mixed. This 
is possible if the rate at which fuel and oxidizer diffuse into each other is 
greater than the rate at which they chemically react with each other. On 
the other hand, a diffusion flame occurs if the diffusion rate of fuel and 
oxidizer is slow compared to their reaction rate. In such a situation, the 
flame exists only in a “reaction zone", where the fuel and oxidizer, which 
are initially separated, mix together through molecular diffusion and 
turbulent dispersion. 

Natural fires are free convection, diffusion flames; their analysis 
involves both solid and gas phase phenomena. After the ignition from an 
external source the combustion in the gas phase feeds heat by radiation and 
conduction to the solid fuel, which in response, generates fuel volatiles, 
which in turn, sustain the combustion. This conversion of a solid phase 
into gaseous volatiles involves the phenomenon of pyrolysis. According to 
Emmons (1980) and Atreya (1983), pyrolysis is one of the three separable 
components in fire phenomena, with other the two being ignition and fire 
spread. Moreover, cellulose maintains its dimensional stability while burning, 
which results in a "classical" diffusion flame enveloping the sample. The 
study of cellulose pyrolysis is necessary in order to understand diffusion 
flames, and the results can be readily applied to diffusion flame theory. 

One of the modern developments in fire science has been the use of 
mathematical models. Such models are used to gain a more fundamental 


understanding of fire phenomena, to interpret standard fire tests, and to 
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design fuel materials or buildings with a prescribed degree of fire safety. 
Diffusion flame models, e.g. Pagni (1980) and Sibulkin et al. (1982), have 
been developed to predict burning behavior such as burning rate and 
extinction limits for a solid fuel with simple geometry such as a vertical 
wall or a stagnation region. These models have also identified the 
material properties which most affect burning. 

Most diffusion flame analyses to date have modeled a simple fuel, 
such as PMMA. For these fuels there is no accumulation of char at the 
fuel surface. Consequently, it can be shown that after an initial transient 
the surface temperature, heat of gasification and burning rate are 
independent of time, as long as the fuel remains thermally thick. On the 
other hand, materials which char, such as wood, show a dependence on time 
related to the growth of the char layer. Sibulkin and Tewari (1985) have 
studied the burning of cellulose, a particularly simple charring material. 

In order to obtain quantitative results from diffusion flame models for 
the combustion of cellulose, a number of materials properties of the 
fuel-oxidizer system must be known. Two of these parameters, which most 
affect burning, are the heat of combustion and the heat of gasification, and 
their ratio is a key factor in burning rate calculation. [trisnticenecn =t OL 
these data that make the study of cellulose pyrolysis very important to 
diffusion flame theories. Pyrolysis and combustion of cellulose have become 


closely coupled topics in fire research. 


ets. 


1.5 The Present Work 


It has been revealed by the literature reviewed in the previous sections 
that_-most | researcher onmaine pyrolysis of cellulosic materials is primarily 
concerned with details of the chemical processes which occur as the materials 
are heated. For this type of work, the appropriate experimental techniques are 
those of thermal analysis such as TGA and DSC, and the typical sample is a 
powdered material weighing on the order of milligrams. Although detailed 
chemical studies are necessary to obtain workable pyrolysis mechanisms, they 
are of little help in actual design of a flame resistant material. 

In the present work, pyrolysis of bulk cellulose in a _ nitrogen 
atmosphere is studied. Sodium hydroxide is used as a fire retardant. The 
materials property of pure and fire retarded cellulose which most directly 
affect its burning behavior, the heat of gasification, is measured, with 
concurrent numerical analysis of one-dimensional pyrolysis. The heat of 
combustion of the volatiles, their stoichiometric ratio and the mass fraction 
of noncombustible volatiles produced in pyrolysis are also obtained as 
experimental results. These data canbe used. as inputs to a~ diffusion 
flame model, and the calculated results can be compared with previous 
cellulose burning measurements. The major objectives of the work are: (1) 
to develop techniques which will be useful in measuring the combustion 
properties of charring materials and which will assist in the development of 
scientifically-based test methods; and (2) to use the experimental results to 
determine the physical mechanisms by which solid phase fire retardants act 
to suppress flaming combustion, which will help to provide a rational basis 


for the design of materials with improved fire safety characteristics. 
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CHAPTER 2 


HEAT OF GASIFICATION OF PURE AND 
FIRE-RETARDED CELLULOSE 


eae | Introduction to Mcasurcments 


AS describedmiineescctionmeico. most rescarch on the pyrolysis of 
cellulosic materials has been primarily concerned with the chemical processes 
which occur as the material is heated, and the typical sample for this type 
of work is a powered material weighing on the order of milligrams. Another 
motivation of the study of pyrolysis is related to its role as one of the 
processes in the burning of solid fuels. The burning of a solid fuel, such 
as cellulose, can be predicted by diffusion flame theory for several simple 
geometries, such as a vertical wall or a stagnation region. In order to 
obtain quantitative results, a number of combustion-related material 
properties must be known. In the present study, a new type of experiment 
is developed, and one of the material properties which most affects burning, 
the heat of gasification, is measured for bulk pure and retarded cellulose 
samples. The measured data will be used as inputs to a diffusion flame 
model. The results of the numerical calculations will then be compared 
with the previous combustion measurements in order to understand how the 
retardants act to suppress fires. 

A heat of gasification is defined as the ratio of the net surface 
heat flux resulting from flaming combustion to the volatile gencration rate 


due to pyrolysis. 


=19a 


In? the he measurements, the basic theory is as follows. The 


definition of the heat of gasification is 


Fl Gnet my 

g aah Tat (2- ) 
m 

where: m = sample mass loss rate per unit area; i.e, the mass flux of 


gaseous products leaving the surface. 


qn net heat flux going into the sample through the front 
surface per unit area and time 
= capes Cares P ahs (2-2) 
where: 1) dab.s = heatfluxabsorbed by sample front surface. 
2) emit = heatflux emitted from sample front surface per 


unit area and time. 


The radiation heat transfer theory for a grey body is used to 


estimate the term: 


" 


emits a €, 9 Che a lies) (2-3) 


where: €, = emissivity of sample front surface = 0.96 


o = Boltzmann constant = 5.67 x 10°8 W/m2 K4 


i 
I 


temperature of sample front surface . 


ambient temperature = 298K 


— 
i 
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3) oe = heat flux due to convection through sample front surface 


per unit area and time. 


The convective heat loss qt could not be measured and therefore, the 


natural convective heat transfer theory for vertical 


planes is followed to 
estimate the term: 


ie = h(T, a hes: (2-4) 


where: h = average value of heat transfer coefficient over sample front 


surface area. 


The effect of "blowing" is approximated by 


- h,&n(1+By) 


Bh 


where B} = blowing parameter 


Cym 
By ape Race 
ho 
: se pe 
hy = average zero mass flux heat transfer coefficient = Ss Nugy 
Nug = average Nusselt number based on sample diameter d 
9/16] 
tes enn 0.492 
Nug = 0.68 + 0.67 Grq Pr 1+ 
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Grq = average Grashof number based on sample diameter d 


S S(T eeacla 


= 
Qa. 
Il 


2 
rive Vr 


The values of Cp» k,, Pr, and v, are taken as constants at T, from 
the physical property table for N,. 
ies Tes 


1.e. T, = reference temperature = Peery ="700K 


‘e) 
i 


p specific heat capacity = 1097 J/kgK 


a 
i] 


thermal conductivity = 0.051 w/mK 


~~ 
4 
I 


Prandtl number = 0.69 


v, = kinematic viscosity = 65 x 10°° m?/s, 


From the equations (2-2), (2-3) and (2-4), the expression for net heat 


flux into the sample front surface is: 


A Po | 4 4 —_ 
Anet = Tabs - £59 (Ts -298 )- h (T, -298) (2-5) 


Thus, the basicy iméasurements. needed -to obtain théeshestsoe 


gasification are the absorbed heat flux, the surface emissivity, the surface 


78 


temperature and the mass loss rate. Detailed descriptions of dap 5, €5. Ts 


- iW é ° 
and m measurements are presented in Sections 2.3, 2.4, 2.5, and 2.6. 
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sen Expcrimental Apparatus 


A schematic side view of the apparatus developed to measure the 
heat of gasification of pure and fire-retarded cellulose is shown in Fig. 2.1, 
while a schematic drawing of the top view of the experimental set-up 
Shown in Fig. 2.2 gives more detailed information about the apparatus. 

The pyrolysis chamber is an extensively modified glass draft shield 
of a top-loading electronic balance. The top of the pyrolysis chamber has an 
added cone through which the gases exhaust. The cellulose sample, in the 
form of a horizontal 2.5 cm diameter cylinder, is held by a sample holder 
made of ceramic foam. The sample and the holder are mounted inside the 
chamber. Because the sample must be enclosed in a chamber having a 
front window, a combination of a quartz plate (acting as an intermediate 
window), radiation shields and cooling coils were added to prevent breakage 
of the front window and to reduce heat flow to the electronic balance. In 
addition, air cooling (not shown) is used to keep the temperature of the 
chamber window and quartz plate below 200°C so that the thermal radiation 
from these surfaces does not affect the measurements. An opaque shutter 
is used to begin the experiment at a well defined time. 

A 2-stage nitrogen supply system is attached to the weighing chamber. 
A valve and tubing arrangement allows the chamber to be purged with pure 
Ny at about 175 &£/min in order to remove air before the pyrolysis 
experiments, and to supply the chamber with pure N> at about 25 2£/min 
during pyrolysis to exclude oxygen from the atmosphere. An inclined tube, 
water manometer is used to monitor the pressure inside the weighing 


chamber to make sure no outside air is leaking into chamber, which shows 
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positive pressure (in the orfer of 1mm H50) during pyrolysis. 

The reason for using a nitrogen atmosphere during pyrolysis is 
related to the desire to simulate the burning of diffusion flames. Both 
experiment (Seshadri and Williams, 1978) and calculations (Sibulkin et al., 
1982) have shown that for many types of diffusion flames the oxygen 
concentration between the flame and the fuel surface is close to zero. Thus 
the gas composition at the fuel surface during burning in air is made up of 
nitrogen and products of combustion (primarily H,0, CO, and CO), plus 
possible traces of oxygen. This atmosphere is simulated by pure nitrogen 
in this experiment. In order to prevent air leaking in, the pyroalveis 
chamber is sealed at the sides by duct tape and platic sheets. The absence 
of oxygen in the chamber was checked periodically by gas chromatograph 
measurements. 

In the experiments, the heat flux from the flame to the fuel surface 
is simulated by a pair of infrared lamps (Research Inc, Model 5305-5A). 
Tungsten filament, quartz envelope lamps giving a_ blackbody spectral 
distribution of energy are used. The lamps are connected electrically im 
series, and supplied with up to a total of 480 volts (240 volts for each 
lamp) through a power controller (Research Inc., Model 646). A control 
panel is set between the power controller and the lamps to control the heat 
flux and the water cooling system. A special mounting holds the lamps at 
an optimum angle (approximately 32° from horizontal) to provide a nearly 
uniform heat flux over the sample front surface. Two heat flux gages 
(Medtherm, 64-20-20T and 64-20-20) are used in absorbed heat flux 


measurement (details are given in Section 2.3). 
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A spectrophotomer (CARY, Model 17) and a Fourier-Transform 
Infrared Spectroscopy (Nicolet, Model 60SX) were used to measure the 
surface emissivity measurement (for details see Section 2.4). An infrared 
pyrometer (Raytek, series SL-300 CC) is used to measure the sample surface 
temperature (for details see Section 2.5). The output of the pyrometer is 
sent to a microcomputer (Apple IIe) for processing. An instrumentation 
amplifier (Interactive Microwave Inc., ADA-AMP) is used for signal 
conditioning of the voltage output of the pyrometer, and then connected to 
the data acquisition interface (Interactive Microwave Inc., ADA-LAB) of the 
computer. Temperature readings are taken every 0.2 seconds (for details 
see Appendix B). 

For the mass loss rate measurements, (for details see Section 2.6), an 
electronic balance (Mettler, PC220) was used. The signal from the balance 
is sent to the computer for data processing through a special interface 
(Microsystems Research Corp. SOFTWEIGH). Mass readings are taken 
every 3 seconds. 

In’ order tO stun thicwexperimcnitsmunder the same conditions, the 
apparatus was checked regularly at specific time intervals. 

The monthly apparatus icluded: (1) heat flux gauge A and B 
(Medtherm, 64-20-20T and 64-20-20) versus multimeter (Fluke, 8840A) at 25 
kW/m? to 40 kW/2 using 5 kW/m? interval; (2) electrical balance (Mettler, 
PC 220) versus standard weight of 50 g; (3) blackbody (Infrared Industries 
Inc., Model 463) with temperature controller (Infrared Industries Inc.. Model 
101B), versus its own thermocouple (Platinum/Platinum 10% Rhodium, Type 


S) at 400°C to 800.C at 50°C interval; (4) ADA-LAB (Apple Ie computer 
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data acquisition interface) reading versus standard voltage source (Analogic, 
Model 3100) at zero to maximum pyrometer voltage output of 4 volts using 
0.5 volt interval; (5) computer (Apple IIe) disk driver speed check. 

The daily apparatus checking included: (1) pyrometer (Raytew 
SL-300) versus blackbody (IR, Model 463) set at 500°C through ADA-LAB to 


both the computer and multimeter (Fluke, 8840 A). (2) the q value 


ab,s 
against the ‘“powér input “to the IR lamps (Research, Inc.) haa 
5305-5A), through multimeter (Fluke, 8840A), by placing het flux gauge A 
(Medtherm, 64-2020T) at the sample location, to make sure the input power 


indicator, RUN POT SETTING, in the control panel corresponds to the 


desired heat flux. 


2.3 Absorbed Heat Flux Mcasurcment 


Heat flux absorbed by the sample front surface is measured using 
the heat flux gauge A (Medtherm, 64-20-20T). The water-cooled radiation 
flux gauge can«be. placed precisely; at the sample location. by,.means of va 
specially made jig. From the slope of the calibration curve supplied by the 
manufacturer, shown in Fig. 2.3., the output of the gauge in mV is 
converted into the heat flux absorbed using the following calibration factor: 
gauge constant = 2.295 W/cm? mV. The heat flux gauge B (Medterm, 
64-20-20) is used as a cross-check of the working gauge. 

A routine check of the heat flux (heat flux gauge A output) versus 
sample height was performed prior to each run using heat flux gauge A 
and multimeter (Fluke, 8840 A). A typical distribution curve of heat flux 


(gauge A output) versus sample height is presented in Fig. 2.4, which shows 
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the sample center absorbs the desired heat flux, which is 40 kW/m2, while the lower 
part of the sample gets more than 40 kW/m2, and the upper part of the sample gets 
less than 40 kW/m. The resultant distribution gives +4% of heat flux variation 
over the sample height. (Previous checks showed the heat flux was uniform in the 


horizontal direction, i.e., parallel tothe filaments of the lamps.) 


2.4 Surface Emissivity Measurement 


According to the radiative heat transfer theory for grey bodies, the 
emissivity, €, iS equal to absorptivity, «a while the absorptivity is equal to one 
minus the reflectivity, p, 1.¢., € =x= 1l- p. Thus, the sample surface emissivity is 
obtained by measuring the reflectivity of the sample surfaces. 

A spectrophotometer (CARY, Model 17) and a Fourier-Transform Infrared 
Spectrometer (Nicolet, Model 60SX) were used to measure the surface relectivity of 
unretarded (Yp = 0) and retarded (Yp = 0.03) samples after different pyrolysis 
times such as 0, 50 and 300s. The spectral ranges of the two instruments are \ = 
0.2 ~ 3.5 um and \ = 2.5 ~ 25 um, respectively, which cover a sufficient range of 
infrared light. The data from the sample surface reflectivity measurements are 
summarized in Table 2.1, from which the average value of 0.04 is obtained. The 


sample surface emissivity is therefore taken to be 0.96. 


9 Surface Tempcraturc Mcasuremcnt 


An infrared pyrometer (Raytek, Series SL 300 CC) is used in the sample 


surface temperature measurement. The pyrometer provides a linear output of 2 to 
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10 volts for a temperature range of 450°C to 1000 °Co. A black body versus 
pyrometer voltage test produced a linear reading of 2 to 10 volts on a multimeter 
(Hewlett Packard, Model #3466A). Thespectral response of the detector is 2.1 to 
2.3 um. Before measurements by the pyrometer are begun, a check is made on its 
reproducibility by focusing it on it on a standard source of black body radiation 
(Infrared Industries Inc., black body radiation source, Model 463). The black 
body itself is checked at 500 °C both by its own thermocouple (for details see 
Section 2.3) and by inserting a thermocouple (Type K, d = 0.003 in.) into its cavity. 
In order to obtain accurate temperature measurement results, a 
calibration curve of the pyrometer voltage output versus sample surface 
temperature measured by a thermocouple is needed. (Direct measurement of 
temperature using thermocouples were not done because attachment of 
the thermocouples interfered with the pyrolysis process.) A series of special 
runs was made for this purpose at qabs = 40 kW/m2, following the 
standard procedure, except that the samples were not painted with black 
heat resistant enamel because of the desire of opCimne a low surface 
temperature. Many difficulties were encountered while trying to attach the 
thermocouple to the sample surface. Several kinds of cements, such as 5 
minute Epoxy, 732 RTV, Sauereisen, and Duco cement, were tried on 
cellulose samples and failed in use because of either thermal or mechanical 
weaknesses. It was decided finally to use Omega CC high temperature 
cement (mix 3 parts filler to 1 part liquid by weight, or 2 parts filler to 1 
part liquid by volume) to glue the thermocouples (Chromel-Alumel, Type K, 
d = 0.003 in.) at the centers of surfaces of unretarded (Yp = 0) and 


retarded (Yp = 0.03) samples respectively. The voltage outputs of both the 
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thermocouple and the corresponding pyrometer are simultaneously recorded 
on a strip chart recorder (Hewlett Packard, Model 7100 B). A schematic 
drawing of a recording is shown in Fig. 2.5. Using a 3rd order polynomial 
curve fit generated on a computer, the calibration curve of sample surface 
temperature versus the pyrometer voltage output shown in Fig. 2.6, is 
obtained. The voltage-temperature polynomial expression corresponding to 
the curve is used in the main computer program to convert the pyrometer 


voltage output into sample surface temperature. 


2.6 Mass Loss Rate Mcasurement 


The mass loss rate of samples during pyrolysis is measured by an 
electronic balance (Mettler, PC 220). The sample is put in the sample 
holder, mounted on the balance and aligned with the pyrometer. The 
weighing chamber is then sealed at its sides using Scotch masker tape. The 
balance is then tared. Once pyrolysis is imitiated, the signal from the 
balance is sent to the computer for processing. In this manner, the mass 


loss rate data are obtained (for details see Section 2.6). 


ed | Results and Discussion 


Heat of gasification measurements were primarily made at an 
absorbed heat flux of qab.s of 40 kW/m. This heat flux gives mass loss 
rates which are comparable to those found by Tewari (1985) for cellulose 
samples burning in air. The typical computer outputs of surface 
temperature, mass loss rate and heat of gasification versus time. for an 
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unretarded sample at the absorbed heat flux of 40 kW/m2 are shown in 
Fig. 2.7 (a), (b) and (c), respectively. It is obvious that surface temperature 
increases rapidly during the first 100 seconds of pyrolysis and then but at a 
slower, Pnate; The surface temperature increases from 700K(427 °C) to 
850K(577 C) during the run. The mass loss rate increases rapidly during 
the first 55 seconds of pyrolysis, reaches the maximum value of 10 x 
10°3kg/m2.s, andwvthen) decreases; afterwards atemancslower) irate, The 
resulting curve for heat of gasification drops rapidly at first and then 
approaches a nearly constant value of hy = 1x 10° J/kg. 

Additional measurements were made at absorbed heat fluxes of 25, 
30; eands5s kW/m? for unretarded samples. The results for different heat 
fluxesusare« shown -in«Fige 2.8,,.Fig.. 2.9, .and, Figs, 2.104 wASmthepheate (lar 
decreases, the values of surface temperature and mass loss rate decrease and 
heat flux decreases, the values of surface tmperature and mass loss rate 
decrease and heat of gasification values increase. While is m" and hg 
kW/m2 to 25 kW/m? makes the biggest quantitive difference for the Teen 
and hy values. Also, the time needed to achieve the maximum mass loss 
rate is affected by the absorbed heat flux; the lower the flux, the longer the time. 

In order to study how a solid phase fire retardant affects the heat 
of gasification hg, runs were made with cellulose samples retarded with 
sodium hydroxide at mass concentrations of 0 to 0.03 (i.e., Yr = 0 to 0.03). 

Fig. 2.11, Fig. 2.12 and Fig. 2.13 show the measured’ ‘values of 
surface temperature, mass loss rate and heat of gasification as functions of 
time for Yp = 0 and Yr = 0.03. Fig. 2.11 indicates that the surface 


temperature for Yp = 0 and Yp = 0.03 increases continuously with time. 


ER = 


The difference between the temperatures decreases with time and reaches a 
constant value after 150 s. Fig. 2.12 shows that for both Yr = 0 and 0.03, 
the mass loss rates increase to peak values and then decline. On the other 
hand, this figure indicates that retardant addition results in a much higher 
mass loss rate value and a faster changing rate of mass loss rate with time. 
The mass loss rate for Yr = 0.03 reaches its peak value at 40 s, while the 
one for Yp = 0 reaches its peak value, which is about 60% of the retarded 
one, at a. later. timesesavours. Also the difference between the two mass 
loss rate values decreases with time rapidly, especially after the peak. Fig. 
2.13. shows that the retardant reduces the value of the heat of gasification 
significantly. While the values of heat of gasification for Yp = 0 and Yp 
= 0.03 both decrease to one-third of their initial values by the end of the 
pyrolysis, the ones for Yp = 0.03 are about 2 times lower than the one for 


YR = 0. The difference between the two h, values decreases with time 


8g 
rapidly, and approaches a nearly constant value after 150 s. As a whole, the 
results show that retardant addition causes increases in surface temperature 
and mass loss rate (Or ‘pyrolysis srate),...and a decrease in the heat’ of 
gasification. It can be seen from the figures that the curves for T, for 
both Yp = 0 and 0.03 have nearly the same shape and are parallel to each 
other except at the very beginning. This is also true for he, but not for m. 
The m" curves show a significant difference in the peak value region; the 
one for Yp = 0.03 is much steeper than the one for Yp = 0. Therefore, the 


increase of mass loss rate plays an important role in the reduction of heat 


of gasification upon retardation. 
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Measured values of surface temperature, mass loss rate and heat of 
gasification as functions of retardant concentration at t = 200 seconds and 
an absorbed heat flux of 40 kW/m2 are shown.in Fig: . 2.14, Fig. 2:)3> ame 
Pca k LG: It can, be seen from Fig. 2.14 that the surface. temperature De 
increases with retardant concentration NOiry On the other hand, AN ~hanges 
with xP in a nonlinear manner. It increases with YR rapidly at first, and 
then keeps increasing but at a much lower rate. In other words, T, becomes 
insensitive to Yp beyond a certain extent. This implies that NaOH is 
acting as a catalyst during pyrolysis, and that only a small amount of 
NaOH is very effective. Fig. 2.15 shows the mass loss rate m" also 
increases with retardant concentration. Similar to the case. for T,, the 
initial small amount of retardant addition indicates a significant influence 


" 


on m". As a result, the heat of gasification, h,, decreases with Yp, as 


g 
shown in Fig. 2.16. Again, a small concentration of retardant causes a large 
decrease in hg. When Yp changes from 0 to 0.03, hy decreases to nearly 
one-third of its initial value. Such a major réduction in the heat of 
gasification with retardant addition is caused by the combination of the 
increase in mass loss rate and the increase in surface heat loss. While the 
increase in T, from 830K at Yp = 0 to 860K at Yr = 0.03 is not large, it 
increases the radiative heat loss from the sample surface enough to have a 
significant effect on aby and hg (for details see Section 2.1). 

The experimental results are contrary to what might have been 
expected, because retardant addition is known to increase the char density 


and presumbly the char conductivity. Heat transfer analysis would 


therefore predict a decrease in surface temperature with retardant addition. 


ache 


This question led to the next phase of the project which was directed at 
measuring noncombustible gaseous products of cellulose pyrolysis (i.e., 


Chapter 4). 
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CHAPTER 3 


COMPARISON OF EXPERIMENTAL AND 
THEORETICAL PYROLYSIS RESULTS 


3.1 Pyrolysis Thcory 


From a conceptual viewpoint, flaming combustion can be divided into 
two. parts: the gas-phase combustion region (the flame) where the fuel 
volatiles are consumed, and the solid phase pyrolysis region where the fuel 
volatiles are generated. The energy released by combustion is transferred to 
the fuel surface as a convective heat flux and a radiant heat flux. The sum of 
the two fluxes, minus a significant heat loss from the hot fuel surface due to 
re-radiation, gives the net heat flux into the fuel. It is this net heat flux 
which causes the fuel to pyrolyze, leading to the production of fuel 
volatiles at a mass flow rate MG w> i.e., the mass loss rate during pyrolysis. 

Three terms are used in the literature to describe the energy changes 
which occur during pyrolysis, they are heat of reaction, heat of pyrolysis 
and heat of gasification. 

The heat of reaction hy @ is defined as the enthalpy of reaction per 


unit mass of volatiles at an ambient temperature T,, for the reaction 
Solid (S) ~ Char (C) + Gaseous volatiles (G) 


so that 


(3-1) 


The definition for the heat of pyrolysis hy is the difference between 
the enthalpy of the products at the surface temperature Ty, and the 


reactants at the ambient temperature T. per unit mass of volatiles so that 
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h. =h +h pe Ne ‘a 
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(3-2) 


Thus hy is equal to hy @ plus the sensible heat needed to raise the 
temperature of the products from T. to T,,. 

For pure cellulose, hy should not depend on the type of sample 
(powered vs. bulk). However, the value of Dy may vary with the heating 
rate, since the products of pyrolysis are heating rate-dependent. The value 
of hy may well depend on the type of sample since for a bulk sample, as 
shown in Fig. 3.1, the initial products of pyrolysis may undergo additional 
chemical reactions as they flow through the hot char layer to the sample 
surface. | 

The definition of hy for use in diffusion flame theory is the ratio 
of the net surface heat flux resulting from flaming combustion to the 
volatiles generation rate due to pyrolysis, that is 


ti Gnet 


g = (3-3) 


m 
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where mG Ww = m". 
An analysis of pyrolysis in a charring material, such as cellulose, 


resulting from an external heat flux (Sibulkin, 1986) shows that the heat of 


La 


gasification is not simply related to the standard heat of reaction during 
pyrolysis. The semi-infinite control volume used in the analysis for a 
charring fuel such as cellulose is shown in Fig. 3.1. Thecoordinate ox) = 
QO represents the fuel surface, and the positive direction of x is taken into 
the fuel. In this analysis, it is assumed that the flame is established at time 
t = 0, when heat flux net begins to go into the fuel. As time increases, a 
As time increases, a pyrolysis wave moves into the fuel converting th original 
solid material (S) into char (C) with the release of gaseous volatiies (G). 

The velocity of the pyrolysis wave, Vp» is related to the mass pyrolysis rate 


Gey by 


ne myles Pc] (3-4) 


An outline of the analysis is as follows: Applying the conservation of 
mass to a differential element of length dx and unit cross-sectional area (see Fig. 


3.1) gives 


" 


G dp 


dx bor 


dm 
(3-5) 


where mc is taken as positive in the negative x direction. 


And applying conservation of energy to the same element gives 


+o |ahe | (3-6) 


Integration of Eq. (3.6) from x = 0 to @ and use of Eq. (3-5) gives 
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After \crelatinguithemiocaltaidensitys p,miand enthalpy h, to the 
corresponding values for the solid and char, the results can be combined 


with Eq. (3-3) to give as a general result for the heat of gasification 


Silico! are 
betes hore _— Iye ae 
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In Eq. (3-8) the first integral on the RHS of the equation gives the 
change in sensible enthalpy between the char and solid while the second 
integral is an unsteady, energy storage term. The analysis thus shows that 
the relationship between the heat of gasification and the heat of reaction 
(or the heat of pyrolysis) is not a simple one. 

In order to gain more understanding of the pyrolysis mechanism and 
as a "check out" of the experiments, a numerical analysis of the pyrolysis of 
pure and fire retarded cellulose was made. A computer program based on 
the method of finite differences has been used to solve Eqs. (3-5) and (3-6) 


to obtain numerical results. The program was developed by Kung (1972) 
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and modified by Tamanini (1976). The complex pyrolysis chemistry is 


replaced by a single global reaction with a rate law given by 


) 
-5, 7 (e-ed A expl-E/RT] (3-9) 


A number of thermal and material properties of cellulose and char 


are needed as the inputs to the program. The results of the calculations 
depend on the values used for the properties. However, accurate values 
of many of these properties are not known. The set of property values 


used in the,calculations are given, in) Table 3.1; ,they are takén) as opese 
estimates from the literature. The thermal conductivities of cellulose and 


char are both taken as the linear functions of local temperature given by 
kWaks- +> k* (T-"Fs) (3-10) 


The values of materials properties, such as the density of char, Pc, the 
thermal conductivity of char at the ambient temperature, KC o and the heat 
of reaction at the ambient temperature, hy ew are, tabulatedsain, Fig. 3:2 64-Fig, 
3.3, and Fig. 3.4, because these properties have proved to be the ones 
which cause the difference between pure (Yp = 0) and retarded (Yp = 0.03) 
cellulose in the calculations. The values of pc came from direct 
measurement, the value of KC @ was taken from the literature, and the hy © 
value came from the DSC measurement results (for details see Appendix C). 

With these inputs, the program can then be used to calculate the 
fuel surface temperature, Te, the mass loss rate, m", and the heat of 


gasification, h,, under different net heat fluxes into the fuel surface, Bac 


s’ 
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(or absorbed heat fluxes dab 5): In order to compare the numerical results 
with experimental data, the heat flux absorbed by the fuel surface was 


taken to be constant at a value of 40 kW/m (see Section .2.7), so that 


q" pos Ga _ q" ; tf 
net ab,s emit,s C,S 


40. ego La) 2) h( TT.) (3-11) 


which is the same as Eq. (2-5) in Section 2.1. As final results, the values 


of Ty m" and hg as functions of time were printed and plotted. 


Bez Theorctical Results 


By using the computer program described in Section 3.1, the 
theoretical results for the pyrolysis of pure and fire retarded cellulose at 
Gab's = 40 kW/m cellulose are obtained. 

Figure 3.2 shows the numerical results of surface temperature Fegaspa 
function of time for pure cellulose (Yp = 0) and retarded cellulose p (Yr = 
0.03). It can be seen that the values of T, for both pure and retarded cellulose 
increase with time, and the T, of retarded cellulose is higher than that of pure 
cellulose, which is qualitatively consistent with the experimental results. 
However, the difference of T, between pure and retarded cellulose is much 
smaller than the experimental one shown in Fig. 2.11. A very important 
phenomenon uncovered by the calculations is that the value of the hcat of 
reaction is different in the two cases: The dominant reactions for pure and 
retarded cellulose are endothermic and exothermic respectively (see Appendix 
C). This result can be used to explain why the value of T, increases with 


sare 


retardant addition (because of the exothermic reaction), . which is 
contradictory to the prediction by heat transfer analysis as mentioned in 
Sechioni2i7: 

Figure 3.3 presents the numerical results of mass loss rate vs. time 
for pure cellulose (Yp = 0) and retarded cellulose (Yp = 0.03). The 
qualitative shape of the. curves, a rapid rise. to .a peak: mass loss rate 
followed by a steady decline, is similar to experimental data (see Fig. 2.12). 
It was confirmed that the addition of the retardant leads to an increase in 
the amass slosseurate A detailed comparison of the theoretical and 
experimental results for the mass loss rate is present in the next section. 

Figure 3.4 gives the numerical results for the heat of gasification as 
a function of time for pure cellulose (Yr = 0) and retarded cellulose (Yp = 
0.03). It shows that addition of the retardant lowers the heat of gasification, 
which is qualitatively in agreement with the experimental data (see Fig. 
2.13). However, the calculated hg values decrease faster than the 
experimental ones, and become constant thereafter. Diffusion flame 
theory predicts that a decrease in h, would promote burning. It was 


8 


expected, therefore, that char-promoting retardants would increase h,. Since 


g: 
the opposite trend was found, further work was carred out (for details see 


Chapter 4). 


aed. Comparison of Experimental and Theorctical Results 


From Fig. 3.2, Fig. 3.3 and Fig. 3.4, it is obvious that the theoretical 


results for both pure and retarded samples agree qualitatively with the 
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experimental results; that is, the surface temperature TT. increases 


S 
with time, m" first increases rapidly to a peak, and then decreases with 
time, and heat of gasification, hg, decreases with time first rapidly and 
then slowly. On the other hand, the addition of retardant leads the 


higher T, and m" values and the lower hy value. 


There are recognizable quantitative differences between the numerical 
and experimental results. The values of surface temperature TS from the 
calculations, shown in Fig. 3.2, are smaller than those from the experimental 
data. The calculated hy values in both pure and retarded cases are 2 to 4 
times higher than the experimental data. A more detailed comparison 
between theoretical and experimental values of the mass loss rate for pure 
(Yp = 0) and retarded (Yp = 0.03) cellulose is shown in Fig. 3.5 and Fig. 
S.C. As mentioned in Section 3.1, the thermal conductivities of cellulose 
and char and the heat of reaction at the ambient temperature, hy eo areethe 
parameters which most affect the pyrolysis calculations, therefore, the 
calculations have were performed for both constant and variable k and hy © 
values. It can be seen from Fig. 3.5 and Fig. 3.6 that for both pure and 
retarded cellulose, the constant values of thermal conductivity and heat of 
reaction lead to results that divergent the most from the experimental ones. 
For the thermal conductivities which vary with the surface temperature 
linearly, in the pure cellulose case shown in Fig. 3.5, the hy value of 
0.5x109 J/kg leads to a peak value which is very close to the experimental 
one, while the Ny value of 0.23x10° J/kg results in a much higher peak 
value but narrows the differences between the "downhill" theoretical and 


experimental data. Similarly, in the retarded cellulose case shown in: Fig, 03.6, 
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the Ay © value of -0.10x10° i/keoieads to the m" values lower than the 
experimental data, while the hy @ value of -0.30x 10° J/kg leads to a higher 
peak value and lower "downhill" values, which indicates a better simulation 
of the experimental results as a whole. 

The Dy @ values used in both comparisons are the best possible 
estimated values available. The comparisons show that for the pyrolysis 
theory used, variable thermal conductivities yield a better simulation of the 
experimental results. 

The numerical calculations confirm the experimental results and show 
that retardant addition increases the surface temperature (this leads to a 
decrease in Gaant and «the j,mass_. loss. ,rate,,..and » decreases, the . heatigor 
gasification. This is contrary to the known fact that retardation increases 
char production, density and presumably the char conductivity; thus there 
should be decreases in both surface temperature and volatile production and 
an increase in heat of gasification. Although it was discovered that the 
value of the heat of reaction is different for pure and retarded cellulose, 
which can be used to explain the increases in T, by retardant addition, 
more information is needed in order to completely understand the 


phenomena. 
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CHAPTER 4 


NONCOMBUSTIBLE GASEOUS PRODUCTS OF 
CELLULOSE PYROLYSIS 


4.1 Introduction to Measurements 


Diffusion flame theory predicts that a decrease in hg would promote 
burning. Therefore, it had been expected that char promoting retardants 
would increase hg. However, our experiments show that retardant addition 
increases the surface temperature and the volatile production rate (i.e., mass 
loss rate), and thus decreases the heat of gasification of cellulose samples. In 
order to explain the contradiction between the theoretical prediction and the 
experimental results, and to explore the mechanism of fire retardation, 
measurements of noncombustible gaseous products of cellulose pyrolysis were 
made. 

Ne analytical gas chromatograph was used to measure the inert 
fraction of the total volatiles. The amounts of noncombustible gaseous 
products, CO, and H 0, in the total volatiles are obtained (the amount of 
CO is also obtained at same time because of the column used). With the 
known inert fraction of the volatiles, the combustible fraction of the total 
volatiles can be calculated, and then used to obtain the heat of combustion 
and stoichiometric ratio of the "fuel", (i.e. the combustible fraction of the 
volatiles). These property values can then be used in combustion 


calculations for burning rates and flammability limits. 


Pals 


42 Experimental Apparatus 


A schematic side-view of the basic apparatus used in noncombustible 
gaseous products measurements is shown in Fig. 4.1. 

In this experiment, -as in the previous heat of gasification 
measurements, a pair of infrared radiant lamps (Research Inc., Model 
5305-5A) are used as the external heat flux source simulating the heat flux 
from the flame to the fuel surface. The same extensively modified glass 
draft of a top-loading electronic balance (Mettler, PC220) is used as the 
weighing chamber, together with the same water-air cooling system and 
2-stage nitrogen supply system. The electronic balance is connected to the 
computer (Apple IIe) to acquire and process the mass data to give the total 
volatiles production rate, 1.e. the mass loss rate, as a function of time. 

An analytical gas chromatograph, i.e. GC, (Carle Instruments, Inc., 
Model 111) is used, together with an integrator (Hewlett Packard, Model 
3380A), to measure the noncombustible gaseous products of cellulose 
pyrolysis. In order to obtain reliable GC results, a special sampling system 
was designed. First of all, based on Heskestad’s study (1981), a system was 
designed to mix the volatile gases which come out from the weighing 
chamber. As shown in Fig. 4.1, the weighing chamber is connected to a 
glass tube where the volatiles are mixed. In order to obtain good mixing, 
an annular metal washer is put at the entrance to the mixing tube. The 
ratio of the washer inner diameter to tube inner diameter is 0.5 as 
recommended by Heskestad for the optimum mixing. 

After mixing, the flow is split in two; i.e, to the GC sampling loop, 


and to a rotameter (Gilmont Instruments, Inc., Flowmeter Catalog No. F1500) 
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used to measure the volumetric flow rate of the volatiles. The rotameter is 
then connected to the main exhaust system through a flexible tube. 

In order to accurately measure the volatile flow rate, a thermocouple 
(Ghromel-Alumel, Type K, ad 290.003" 1m) “is ‘inserted into “the ‘mixing ‘tube. 
The thermocouple is connected to a digital temperature indicator (Omega 
Engineering, Inc., Model 199) which reads the temperature of the volatiles 
passing through the rotameter. These temperature data are then used in the 
calculation of the volatiles flow rate. (The flow rate cannot be read directly 
from the calibration chart, because the volatiles temperature, ranging from 
56°C to 72°C is far above the standard temperature, i.e. ZINE, at which the 
calibration chart is obtained.) 

A vacuum pump (Cole-Parmer Instrument Co., Air Cadet, Model No. 
TII0°40)" "IS used’ to “create “the “presstire “difference needed “to draw’ the 
volatile gases through the GC sampling loop. (Before injecting volatiles 
sample into GC columns, the pump is turned off to let the pressure inside 
the GC sampling loop recover.) The pressure in the GC sampling loop is 
monitored by a mercury manometer. The volatiles flow rate through the 
sampling loop is measured by a rotameter (Emerson Electric Co., Sho-Rate 
150, Model 1355). A pump bypass is used to adjust the flow rate through 
the sampling loop to a value of 180 m2/min during GC sampling. The 
time needed for volatiles to move from pyrolyzing sample surface to GC 
sampling valve has been estimated so that during the experiment, the 
volatiles inside the sample loop can be injected into the GC column at the 
proper time, which is directly related to the time used for the mass loss 


rate measurements. ‘The GC used helium as the carrier gas. The helium 
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flow was 22 m&/min at a pressure of 42 psig. Before the measurements, 
the helium flow was run through the GC columns by using the valve switch 
to get air out of the sampling loop (because of a possible overnight 
leakage). The column temperature was set at 60°C. A schematic drawing of 
GC valves and columns is shown in Fig. 4.2. There are 3 columns in the 
GC, column 1 (6 -80/100 Porapak Q) is used to measure the noncombustible 
volatiles, CO, and H,0, while column 2 (5-80/100 Molecular Sieve 5 Auoag 
used to measure O54, N>5, and CO and column 3 is used as the reference 
column. 

The amount of volatiles injected into the GC columns is | mf&. 
After injection, the positions of the GC valves are controlled by a timer. 
The volatile sample is sent either through column 1, or column 2 or column 
1 plus column 2. The timing sequence has been chosen so that CO, and 
H50 are separated by column 1, and O>, N» and CO are separated by 
column 2. A TCD (thermal conductivity detector) was used. The output 
of. the GC is sent ie the integrator which obtains the area of the sample 
peak and calculates the area percentage of each peak. The total time of 
analysis was 30 minutes per sample. The GC was calibrated using a bottle 
of calibration gas and a mixture of laboratory air and N». Details of the 
GC calibration are presented in Appendix D. 

In this experiment, as in the previous one, some of the apparatus 
were checked regularly in order to keep the experimental conditions 
consistent. 

The monthly apparatus checks included: (1) heat flux gauge A 


(Medtherm, 64-20-20T) versus heat flux gauge B (Medtherm, 64-20-20); (2) 
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electronic balance (Mettler, PC 220) versus standard weight of 50 grams; and 
(3) computer (Apple IIe) disk driver speed checking by using the Master 
disk. 

The daily apparatus check was made of the IR lamps against the 
heat flux gauge A by using a multimeter before and after each run. The 
purpose of (this wasmtow insuremanarmince UN POT SELTING. in the control 
Panel corresponds. to minemmucsitcumicatriux.. and that the heat flux is 
uniformly distributed over the entire sample surface (for details see Section 
geal 

Since the total GC measuring time is 30 minutes per sample while 
the total cellulose pyrolysis time is 5 minutes, only one GC sample can be 
faken. during. cach irin: Based on the hy measurement results, it was 
decided that the volatile products of cellulose pyrolysis at 90, 180 and 270 
seconds will be measured by the GC. (The time needed for volatiles to 
move from the pyrolyzing sample surface to the GC sampling valve, 


estimated as 15 seconds, was taken into account.) 


4.3 Results and Discussion 


A series of runs were made to measure the noncombustible gaseous 
products of pyrolysis for pure (Yp = 0) and fire retarded (YR = 0.03) 
cellulose samples. The volatiles coming out from the sample surfaces were 
sampled by the GC at 90, 180 and 270 seconds. The GC results were 


recorded by the integrator. 
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For pure cellulose samples, the GC analysis showed peaks for CO», 
H,0, Oy and N>, while for fire-retarded samples, the GC results included 
CO. An example of the GC output from the integrator is shown in Fig. 
4.3 for a fire retarded (Yp = 0.03) sample at 180 s. It can be seen that the 
integrator starts recording at the moment the volatile sample is injected into 
the GC columns. It is known from the GC calibration run that the peaks 
recorded af '2.1945.820014:05,0)4:74.and 20.72eminutesigreprésent CO>,> HAG; 
O5, N> and CO, respectively. The other small peaks are just the responses 
to the GC valve switches from clockwise to counterclockwise. The GC 
analysis is completed at 30 minutes as preset by the GC operating program. 
After the recording is stopped, the GC results are tabulated to give the 
recording time (RT), area (in counts) and area % (in count %), where the 
column "Type" explains the characteristics of the individual peaks . The 
"T" means a peak with "tangent skimmed base", and the "M" means a peak 
which "merged with the preceding peak". It can be seen. from the data 
that among the total volatiles, the area % of CO5, HO and CO are 0.21%, 
0.44% and 0.14%, respectively. Also, there is 99% N» in the volatiles, which 
agrees with the estimated value. The peak representing N» is chopped off, 
because it has exceeded the recording limit. The amount of O>5, 0.02% is 
negligible, which shows that within the acceptable experimental error, the 
cellulose sample surface can be considered surrounded by nitrogen, as in the 
"real" diffusion flame situation. The density of the total volatiles can be 
approximated by the density of nitrogen. 

As described in the previous section, the masses of CO>4, HzO and 


CO are calculated from the GC analysis results, which then is divided by 
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the total volatiles production rate which is measured simultaneously, leading 
to the results for mass fractions of CO? H5,0O and CO. The? {iilar resuits 
of the measurements of noncombustible gaseous products of cellulose 
pyrolysis are summarized in two figures. 

The mass fractions of CO, H,0 and CO vs. time for pure and 
fire retarded samples at canes 40 kW/m2 are shown in Fig. 4.4. In this 
figure, the values of Yco2 YH20 and Yco, with the error bars, are 
plotted against the pyrolysis time for both the Yp = 0 and Yp = 0.03 cases. 
The values of Yco 7 increase with time at about the same rate for Yp=0 
and 0.03. For Yy1709, the values for Yp = 0 increase with time at a faster 
rate than those for Yp = 0.03. The values of Yco for Yp = 0.03 
increase very slowly with time, while there is no measurable amount of CO 
for Yp = 0. The most significant results of the measurements are the 
effects of the retardant on the amounts of CO, and H,0 produced. The 
figure shows that the values of Yco > for Yp = 0.03 are about 4 times 
higher than those for Yp = 0. Also the values of Yyyx09 for Yp = 0.03 
are about twice those for Yp = 0. 

The mass fractions for CO,, H,O and CO vs. retardant 
concentration at dab.s = 40 kW/m? are shown in Fig. 4.5. In this figure, 
the values of Yco7, Yy29 4nd Yco with error bars are plotted against 
retardant concentration for a pyrolyzing time of 180 seconds. Fig. 5.6, 
shows that the values of Yco2 YH20 22d Yco increase with the retardant 
concentration. It also shows that Yco?2, YH20 and Yco all increase very 
rapidly as Yp changes from 0 to 0.01. Over this Yp range, Yco 


nearly increases to its maximum value and remains almost constant 
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thereafter. The value of Ycg increases to 80% of its value for Yp = 0.03 
and then increases at a much slower rate. On the other hand, Yq0 
increases to its maximum value at about Yp = 0.015, and then starts to 
decrease. These results imply that the most significant changes in Yeo), 
Yy290 and Yco are caused by a small amount of retardant addition. 


Therefore the retardant NaOH is speculated to be acting as a catalyst. 


Both figures show that the Yyyx9 values have larger errors than the 
Yco2 and Yco values. From the Fig. 4.3, it can be seen that the peak of 
the H5,O curve has a long tail which is significantly different from the 
other peaks. Therefore, cutting off the H50 peak tail, which is executed 
by the integrator automatically during the GC analysis, can be considered as 
a possible source of error in the HO measurements. 

The results for Yeo and Yyy7o9 are very valuable in explaining 
the puzzle left by the hy measurements, because even though the addition of 
retardants decreases the he values unexpectedly, it also increases. the 
- noncombustible gaseous products, CO, and H,0O, significantly during the 
cellulose pyrolysis. The physical mechanism by which the fire retardant 
NaOH acts to suppress flaming combustion has to be explained by the 


combination of the two aspects. 
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CHAPTER 5 


PRODUCTS OF CELLULOSE PYROLYSIS 


5.1 Introduction to Analysis 


Diffusion flame theory shows that the fuel-related material 
properties which must be determined to predict burning behavior are: (1) he, 
the’ Néat Of Weasification:. 16. athe neat input per “unit mass “of volatiles 
produced; (2) Ypy, the mass fraction of combustible volatiles in the total 
volatile products of pyrolysis; (3) v,, the stoichiometric oxygen/fuel ratio of 
the combustible volatiles; (4) Ao cv the heat of combustion per unit mass of 
combustible volatiles. Of these four properties, hg has been measured, as 


and h can 


described in Chapter 3. The other three quantities, Ypry, V, rey 


be determined as follows from using the experimental results. 
The mass fraction of combustible volatiles in total volatiles, Yp7, is 


defined as: 


™(CO, + H50) 
\ Fl) 3 ihre e 6. CO>) 7. S50 a 
Myo] 


therefore, once the values of Yco2 and Yq 0 are calculated from the GC 
analysis results, the values of Ypy are known. In order to obtain the 
stoichiometric oxygen/fuel ratio of the combustible volatiles, v., the C, H, O 
composition of the combustible volatiles are calculated first from the 


following: 
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where: Yicy = #™mass fraction of species "1 6) inteothedicombustiole 
volatiles. The, subscript 31. representswaG. Hie 


respectively in this section. 


m; cy = mass of species i in the combustible volatiles per 
unit time. 

nee = mass of the combustible volatiles per unit time. 

Y; vol = Mass fraction of species i in the total volatiles. 


mass fraction of species i in CO>. 


ig 
Qa 
S 
ul 


Yco, = mass fraction of CO, in the total volatiles. 
at HO = mass fraction of species i in H50. 
YH,O = mass fraction of H5O in the total volatiles. 
Yer = mass fraction of combustible volatiles in the total 
volatiles. 


The values of Yizco2 and Yi, H20 are obtained by using’ basic 
chemistry. The value of Yco2, YH20 29d Yer are calculated from the 


GC analysis results. 
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The values of Yi vol are calculated by using the following formula: 


My vol 
Yi vol oes : 
Myo] 
Yivcell ~ Yi,char Ychar 7 
= (=>) 
yi Y char 
where m;,,; = mass of species i in’ the total volatiles per unit 
time. 
nisi = mass of total volatiles per unit time. 
Meoeell = mass fraction of species iin the cellulose samples. 
Yi char = mass fraction of species 1 in char. 
Yohar = mass fraction of char in the cellulose samples. 


The values of Yicell are  cakculated’ for** both’ pure “and retarded 
cellulose samples, using the molecular formula C6H 10995 for pure cellulose 
samples. The values of Yicnat are measured for both pure and retarded 
samples by using a elemental analyzer (Perkin-Elmer, 240 B); details about 
the measurement are given by Frendi (1986). TREY Char avalues#ror "pure 
and retarded samples are obtained by making a special series of 
measurements. The cellulose samples are cut at the front end to make 


disks 2-3 mm (thick@and @iesinie(2.54 cm) in diameter, These disks are 


then attached to cellulose samples by a specially made holder. The cellulose 
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samples with disks are pyrolyzed following the same procedure as in he 
measurements, except using an increased heating time in order to let the 
Section 3.3, except using an increased heating time in order to let the disks 


become completely charred. The values of Y are calculated by using 


char 
the mass data of the disks before and after charring. 
Once the C, H, O composition of the combustible volatiles becomes 


known, the values of v,. are calculated by using the following stoichiometric 


equation for complete combustion: 


CxHyOz + 102 = x CO) + 5 H30 


where C,H,O, is the combustible compound. 


From the balance of the equation: 


] y y Z 
nm = ——(|2xX + —-Zz| =x + —- — 
2 2 4 ¢) 
For a unit mass of Cx HyO,, 
YCoov “i YO,cv 
1 1V5 Se ahaa iy 
so that 
YCcv YH ov YO, cv 
feet eh 5 a 
12 4 Fe 
and 
Vi mse Dp (5-4) 
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The values of the heat of combustion per unit mass of combustible 


volatiles he oy can be calculated from 
Ae vol 
Aoey ir (D-9) 
YET 
where he, vo|> the heat of combustion per unit mass of total volatiles, has 


been measured by A. Frendi (1986) using a semi-microcalorimeter (Parr 
1421), for pure and retarded cellulose samples ranging from Yr=O to 0.03. 

The values for hg, Morcys Vg nd Ne ey obtained. aLas 
described above, are used as inputs to a computer program to predict 
the burning behavior of cellulose. 

It is interesting to look at the decomposition of cellulose into char, 
inert gases and combustible volatiles. The amount of char per unit mass of 
cellulose Yonar is from the measurements, as mentioned earlier. The 
amount of inert gases per unit mass of cellulose Yinert 2nd the amount of 


combustible volatiles per unit mass of cellulose Ycomb,vol are calculated 


from the following formulas: 


. Minert 
inert ~ 
M cell 
ee (1 “ errr) (1 = eha re (5-6) 
and 
Ycomb,vol = 1- Yohar > Yinert (5-7) 
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The results are tabulated in Table 5.7 andsplotted-inwkig one 


52 Results and Discussion 


During cellulose pyrolysis, a certain amount of char is produced. 


The value$”" Of "mass traction. o1. char “in” the-*cellugose  sanipices.ae are 


char? 
obtained by charring the thin disks for pure and retarded cellulose samples, 
under the same conditions as described in Section 2.2. The results of these 
measurements are tabulated in Table 5.1. It is obvious from the table that 
the production of char from cellulose pyrolysis increases as the retardant 
concentration increases. However, the mass concentration of char in the 
cellulose samples, Yop4,, does not increase linearly with the retardant 
concentration Yp. As Yp increases, Yop4, imcreases rapidly and then slows 
down approaching a maximum value. 

Table 5.2. gives the elemental composition values for cellulose 
according to the formula CgH QOs. It can be séen that C and O°  dieuene 
aaiot species in cellulose, each composing about 40% of the cellulose by 
mass. 

The elemental composition of char are given in Table 5.3. The 
analysis are described in the last section. It was assumed that all Na 
remains in the char, the data for C and H were obtained from direct 
measurements, and the value for O was found by subtraction. The table 
shows that the mass fractions of C, H and O vary only a little as the 


retardant concentration increases from 0 to 0.03, and can be considered as 


nearly constant. As Yp increases from 0 to 0.03, the values of mass fraction 
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of C, H and O composition of char remain around 80%, 4% and 16%, 
respectively. Therefore the data in Table 6.3 lead to the conclusion that the 
retardant addition does not cause significant effects on the C, H and O 
composition in char. 

The elemental compositions of the volatiles in Table 5.4. are 
obtained by using formula (5-3) shown in Section 5.1 and the data in Tables 
oiggl A cade op A § oD ak Tavlicwesecnows tat ine Values of the C, H and. O 
composition of total volatiles change with the retardant concentration. As 
+62) increases from 0 to 0.03, the mass fraction of C in the total volatiles 
decreases from 41% to 29%, an approximate decrease of 30% of the initial 
value, which is presumably because more of the original C ends up in the 
char, as shown in Table 5.1. On the other hand, the mass fraction of H 
and O increase with increasing Yp. The total increases of the H and O 
in the total volatiles amount to 10% and 20% of the initial values, 
respectively. The fact that the C, H and O composition of the total volatiles 
changes with the increasing retardant concentration shows the Chee Ctsms. 
retardant addition on the elemental composition of total volatiles. 

In order to obtain a better understanding of the physical mechanism 
by which the fire retardant NaOH acts to suppress flaming combustion, data 
on the noncombustible products of cellulose pyrolysis were obtained using 
the GC results. The details about the calculations of the mass fraction of 
CO, and H5O in the total volatiles, Yoo 7 and Yy 0, are presented in 
Appendix D and shown in Figs. 4.4 and 4.5. The values of the mass 
fraction of combustible volatile in total volatiles Ypy are calculated from 
formula (5-1) wT Abe measured values of Yoo? and Yq 0. Tapiew>.o 
gives the data on the noncombustible gaseous products of cellulose pyrolysis. 


ess. 


It shows that as Yp increased from 0 to 0.03, Yeo increased 3 times 
(from 0.06 to 0.24), and Yyyo9Q nearly doubled (from 0.25 to 0.41). As a 
result, Ypqz decreased to half .of, -its initial, value, (from..0.69  tows0 sae 
Therefore, Table 5.5 clearly indicates that redardant addition causes. Ja 
decrease of gas phase “fuel” in volatiles. 

Table 5.6 shows the elemental composition and the stoichiometric 
oxygen/fuel ratio of the combustible volatiles. The values of the mass 
fraction of C, H and O in the combustible volatiles are calculated by using 
formula (5-2) and the data in Tables 5.4 and 5.5, while the values of 
stoichiometric ratios of the combustible volatiles are calculated by using 
formula (5-4), as described in the previous section. It_ can ,be, seen stom 
the table that the C, H and O composition of the combustible volatiles 
changes with increasing retardant concentration in a different way from 
that for the total volatiles. Thesvalue cot YCov first increases as Yp 
increases, reaches a maximum value at Yp=0.015 which is 30% higher than 
the minimum one corresponding to YRp=0, then decreases to a value which is 
about 10% higher than the minimum one. The value of YC oy, om the other 
hand, continuously increases, from the initial value for Yp=0, to a 
maximum value for Yp=0.03, which is 40% higher than the initial one. 
While the Yo oy value first decreases as Yp increases, reaches a4 minimum 
value, which is just half of the initial value, at Yp=0.015, then increases to 
the value, which is about 80% of the initial one, at Yp=0.03. The effects of 


these three factors result in the v, data. Similarly to Yooy, the v 


5 value 


also first increases as YR increases, reaches a maximum value, which is 


Se 


about 44% higher than the initial one at YRp=0.015, then decreases to a 
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value which is 25% higher than the initial one at, Y p=0.03. The very 
important fact here is that the value of Vs increases with increasing ay 
which implies that in an environment with certain concentrations of oxygen 
and combustible volatiles, the addition of retardant NaOH simply makes 
cellulose harder to burn. 

As a summary, in order to give a clearer picture of the results of 
cellulose pyrolysis, the data on mass fraction of the major products of 
cellulose pyrolysis per unit mass of cellulose are tabulated in Table 5.7 and 
plotted against retardant concentration in Fig. 5.1. maple 5-259nd .Ficuso.! 
can be used to explain how the retardant NaOH acts to suppress burning 
through two important factors. First, as retardant concentration increases, 
both the mass fraction of char and inert gases increase causing the mass 
fraction of combustible volatiles to decrease. Therefore, one can conclude 
that the retardant addition reduces the amount of the "fuel" in the volatiles, 
and, thereby, suppresses burning. Second, Fig. 5.1 shows that 90% of the 
mass fraction changes occurs as Yp increases from 0 to 0.015. This 
suggests that retardant NaOH is acting as a catalyst to speed up and 
increase the formations of char, CO, and H,0O at the expense of 
combustible volatiles, as discussed in Section 1.1. Therefore, it appears that 
only a small amount of NaOH is needed for fire retardation of cellulose. 

The overall effects of NaOH addition on cellulose pyrolysis and 


combustion are summarized in the next chapter. 
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CHAPTER 6 


SUMMARY OF RETARDANT EFFECTS 
ON BURNING OF CELLULOSE 


The motivation for this study of cellulose pyrolysis is related to its 
role as one of the processes in the burning of cellulose. In particular, the 
present study is directed toward determining the effects of the solid phase 
fire retardant, NaOH, on the burning of cellulose. The, sresuitew ane 
summarized as follows: 

(1) °The® theory “of diffusion ‘flames’ "Shows Wthdte am viriporiems 
parameter relating gas-phase combustion and solid-phase pyrolysis is the heat 
of gasification hg, which is the heat input per unit mass of volatiles. A 
new technique for the study of bulk cellulose pyrolysis has been developed 
to measure surface temperature, pyrolysis rate, (mass loss rate), and heat of 
_ gasification... [The ¢xperiméntal results indicate .that> addition.) of a4 tire 
retardant (NaOH) to cellulose increases the surface temperature, which in 
turn decreases the net heat input, and also increases the pyrolysis rate, 
(while some previous researchers reported decreased pyrolysis rates after 
adding retardants); these effects result in a lower heat of gasification. 
Analysis of pyrolysis wave propagation in a char-forming material such as 
cellulose shows that there is no simple relationship between the heat of 
gasification and the heat of reaction. Numerical calculations for a pyrolysis 
wave propagating into cellulose using reasonable estimates of the material 
propertics were performed to generate theoretical predictions for surface 


temperature, pyrolysis rate and heat of gasification. The experimental data 
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are in qualitative agreement with the theoretical results, but there is a 
significant difference in the magnitudes of the apparent heat of 
gasification. 

(2) Diffusion flame theory also shows that the mass fraction of 
combustible volatiles in the total volatile products of pyrolysis, Mepwithe 


Stoichiometric oxygen/fuel ratio of the combustible volatiles, v., and the 


5° 
heat of combustion per unit mass of combustible volatiles, Ne evs are) the 
fuel-related material properties which must be determined to predict burning 
behavior. Based on the assumption that the only noncombustible gases in the 
volatile products of pyrolysis are CO, and H»50, a GC was used, (together 
with the experimental apparatus for the he measurements) to determine the 
amount of CO, and H>0O in the volatiles. Calculations using the GC results 
show that the addition of NaOH lowers the Ypry values. The data from 
the GC measurements and char analysis also show that the NaOH addition 
causes an increase in the v, values. Although Frendi (1986) has shown 


that the heat of combustion of the total volatiles, h decreases with 


c,vol 
increasing NaOH concentration, the combined effects of Ae vol and Yrr 
lead to increased Ne ey values with NaOH addition. 

(3) The previous experiments by Tewari (1985) showed that the 
retardant NaOH acts to suppress the burning of cellulose. In the present 
work, the retardant NaOH is found to have a dual effect on cellulose 
burning. Table 6.1 shows the material properties used in diffusion flame 
calculations. NaOH acts to retard (or suppress) burning by: (1) increasing 


surface temperature, thereby decreasing the net heat input to the solid fuel 


surface; (2) decreasing the mass fraction of combustible volatiles in the total 
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volatiles, thereby reducing the gas phase "fuel" during burning; and (3) 
raising the stoichiometric oxygen/fuel ratio of the combustible volatiles, 
thereby making cellulose harder to burn in certain oxygen-containing 
environments. On the other hand, NaOH acts to promote burning by: (1) 
lowering the heat of gasification, thereby increasing the production rate of 
combustible volatiles for a given amount of net heat input to the solid fuel 
surface; and (2) raising the heat of combustion per unit mass of combustible 
volatiles, thereby increasing the heat released per unit mass of gas phase 
epuel 

In order to determine the complete effects of NaOH on the burning 
of cellulose and the detailed physical mechanisms by which solid phase fire 
retardants act to suppress flaming combustion, further studies regarding the 


physical and thermal properties of the char are needed. 
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APPENDIX A 


SAMPLE PREPARATION 


In order to obtain reproducible results for the pyrolysis of bulk 
samples of cellulose, first of all the samples have to be prepared by 
following a standard procedure precisely so that the properties will not very 
from one sample to another. On the other hand, a sample has to have the 


following principal characteristics: 


1. A homogeneous structure to assure homogeneous properties. 

2. A smooth and crack-free surface. 

3. Sufficient mechanical strength to survive experimental 
Operations without breaking. 

4. For a retarded sample, the retardant should be distributed 


uniformly throughout the sample. 


Based on Tewari’s (1985) studies, it was decided that for a sample 
made of a total of 20 g of cellulose, a mix of 8g (2/5 by weight) of filter 
paper pulp (Schleicher and Schnell, No. 289, 0.01% ash content) and l2g 
(3/5 by weight) of cellulose powder, (Whatman, Type CF11, 0.015% ash 
content) is the most suitable combination in terms of surface smoothness and 
mechanical strength of the samples. For retarded samples, the amount of 
retardant, NaOH, is calculated by multipylying Ya, (retardant concentration 
in the cellulose slurry) by 20 g (total amount of cellulose), where Ya can 
be obtained from the calibration curve in Fig. A.] through Yp (retardant 


concentration in the sample). It was also decided that the shape of the 


Gt. 


sample should be a cylinder, 2 in. (5.08 cm) long and 1 in. (2.54 cm) in 
diameter, with flat surfaces on both ends. 

The apparatus for making the cellulose samples is shown in Fig. A.2. It 
consists of a plexiglass cylinder, 24 in. (61 cm) in length and 1! in (2.54 cm) 
in internal diameter, with a piston at the top. The piston can be pushed 
down by a pneumatic press at a controlled speed and pressure. There are 
two porous metal filters on the piston bottom and cylinder bottom, 
respectively. The cylinder is placed in a stainless steel cup. 

A total of 300 ml of distilled water (including the wash water) is 


used in making each sample. The following procedure is adopted: 


1) Weigh out 12 g of cellulose powder and 8 g of cellulose pulp on an 
electronic balance (Mettler, PR700). Put the two kinds of cellulose in a 500 
ml beaker, add 200 mi (175 ml for retarded samples) of distilled water from 
a measuring flask, and stir the mixture with a stirrer. For retarded 
samples, weigh out the desired amount of NaOH on an electronic balance 
(Mettler, H80), place it in a 100 ml beaker, add 25 ml of distilled water, 


and dissolve the retardant by Stirring. 


2) Put 100 ml of distilled water in a wash bottle to wash the 
container surfaces which have been in contact with cellulose (or cellulose 
and NaOH) slurry and NaOH solution in order to collect all the original 


materials. 


3) Pour cellulose slurry (or cellulose slurry and NaOH solution) into a 
blender, wash the beaker containing NaOH solution, and run the blender for 


3 minutes to obtain a well-mixed cellulose (or cellulose and NaOH) slurry. 
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4) Pour the contents of the blender into the 500 ml beaker, and wash 
the sides of the blender. Then pour the cellulose (or cellulose and NaOH) 
Slurry into the cylinder of the sample making apparatus, as shown in Fig. 
2.1., wash the beaker containing the slurry and pour the remaining wash 


bottle water into the cylinder. 


5) Adjust the control valves on the pneumatic press system to let the 
piston be pressed down with a pressure of 110 psig (8 kg/cm?) and an 
average speed of | in/min (2.54 cm/min) to squeeze the water out into the 


cup. 


6) Stop the piston at a calibrated height to produce a 3 in. (7.6 cm) 


long sample. 


7) Drain ail (the water fromethe® cylinder’ by inverting it, take the 


sample out of the cylinder with the help of a special jig. 


8) Dry the sample in a laboratory oven (Blue M, SwW-11TA) for 5 


hours at 85°C to remove moisture. 


9) Cut the dried sample to 2 in. (5.1 cm) length and measure the 


density. 


10) Spray the sample front surface with a black hcat resistant enamel, 
(If the cellulose samples were pyrolyzed in their natural condition, the 
surface emissivity would change with time until the surface was fully 
charred. In order to have a known value of surface emissivity, the samples 


surface is sprayed with the high temperature black paint). 


Wan 


11) Dry the sample in the laboratory oven for ! hour at 105°C to 
remove the final traces of moisture, and then keep it in a plastic dessicator 


(Bel-Art) filled with dessiccant. The sample is then ready for use. 


The samples thus produced have a density of 0.55 + 0.01 g/cem?, 
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APPENDIX B 


DATA ACQUISITION AND PROCESSING 


B.1 Temperature Signal Conditioning 


The voltage output from the pyrometer inputted into the ADA-AMP 
signal conditioner is digitized in an A/D convertor, which gives a running 
digital count of the input voltage. The following are the relations between 


the input voltage (V) and the digital count (DC). 


FS 2(GN+1 
Lory ai iM +1)] ra 


(SW GN)(OP-AMP GN) 


which gives 


_ 2047 _V(SW_GN)(OP- AMP GN) 
FS[2(GN+1)] 


where: 


V = Actual input voltage in volts; maximum +10 volts; 
+ 4 volts in this case. 
DC = Digital count of the input voltage; equals +2047 


at the peak of the full scale voltage range (+FS) 


in the ADA-LAB. 
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FS = Jumper selectable voltage range in the ADA-LAB; 
can be +0.5 v, +1 v, +2 v and +4 v. +4 volts in 


this case. 


GN = Programmable attenuation index in ADA-AMP; range 


fromeo,to.9, 0 in this case: 


SW GN 


Switch selectable gain in the ADA-AMP; can be 0.1, 


1, 10, and 100; 1 in this case. 


OP-AMP GN 


Operational amplifier gain in the ADA-AMP; 


continuously variable from ! to 10; 2 in this case. 


If the input voltage range is known approximately, a judicious 
selection of the above parameters can be made which allows the A/D 
convertor to use its full range, that is, DC as close to +2047 as possible. For 
example, in the pyrolysis experiments, FS, GN, SW GN and OP-AMP GN are 
set as +4 v, 0, 1 and 2, respectively, for the extreme case in which V could 


be the limited value of +4 volts, the equation (A-2) gives: 


os Wage A ae 
+4 x 2 


DC 2047 


2047 


It shows that the gain settings used lead to the maximum utilization of the 


A/D convertor. 
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B2 Tempcrature and Mass Sampling 


Fig. A.1 illustrates the way temperature and mass are sampled. 
It is essential to leave at least a one second margin between the previous 
mass sampling point (say 3 seconds) and the start of the temperature 
sampling point for the next mass sampling point (say 6 seconds). atts 
is because whenever the SOFTWEIGH interface is activated to sample 
mass, the computer and hence all the peripheral interfaces (including 
Ouick -1/0), with the? exceptionmotmstne, Quick 1/0 timer, are. disabled for 
about 1 second. This one second is the time required by the balance to 
digitize the analog signal into an 8 bit parallel signal, and then convert 
it into a serial signal to transmit it through the balance cable. Finally, 
the serial signal is converted again into an 8 bit parallel signal which 
then is processed by the SOFTWEIGH interface. 

Although the computer code can handle arbitrary values of mass 
and temperature sampling rates, because of the reason above, there is a 
limited relationship between the two sampling rates, i.e, a 1 second time 
gap between the mass sampling point and the next temperature sampling 
point. Also, because of BASIC Language Software limitations, a sampling 
interval less than 0.2 seconds is not possible. Therefore, in minimum 
mass sampling time allowed is 3 seconds if 10 temperature samples are 
taken by using the minimum 0.2 seconds sampling interval. A shorter 
mass sampling time, eg. 2 seconds is possible if only 5 temperature 


samples are taken at 0.2 seconds sampling interval. 
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B.3 Tempcrature and Mass Data Processing 

Because of fluctuations in . the -mass readings .as sent to ine 
computer, it is not possible to calculate the mass-loss rate using the oraw 
data. It was decided to run a least squares curve fit through the raw 
mass data, and calculate a curve fit to get the mass and mass loss rate 
at different times. The same strategy was adopted to smooth the raw 
temperature data obtained from the pyrometer, which also. exhibits 
significant fluctuation. 

To generate input data for the mass regression scheme (or curve 
fitting), mass values from the balance, along with corresponding time data, 
at 19 points (starting from 3 seconds) are taken. After fitting a second 
order curve through these points, a "smoothed" mass value is calculated at 
the midpoint of the data, and the corresponding mass loss rate is obtained 
by differentiating the curve fit at that point. Next, again mass and time 
data at 19 points, starting from 6 seconds this time, are chosen, a second 
. order curve fit is done, and mass loss rates are calculated at the curve 
midpoint, i.e. at (6+60)/2 = 33 seconds. This. process is. continued and 
"smoothed" mass and mass loss rate values are obtained every three seconds, 
starting at 30 seconds, and ending at 270 seconds for a run of 300 seconds. 
The curves of mass and mass loss rate versus time are thus obtained. 

The value of the temperature at each mass calculation point (30, 33, 
36 seconds, etc.) is obtained by fitting a second order curve through 10 
temperature raw data readings from the pyrometer. That sis. in order te 
obtain a "smoothed" temperature value at 30 seconds, temperature readings 


from the pyrometer are collected at 0.2 seconds interval, starting at 28 
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seconds, and) a second) ordémecuive 1s Litted through these 110 points. The 
temperature at 30 seconds is thus obtained (see Fig. B.1). This process is 
continued and a smooth curve of sample surface temperature versus time is 


obtained. 
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APPENDIX C 


DSC MEASUREMENTS OF HEAT OF REACTION 


Heat of reaction measurements were made for unretarded (Yr = 0 and 
retarded (Yp = 0.03) samples. A differential scanning calorimeter (DSC) 
was used (Dupont 9900 DSC/DTA system). The DSC was calibrated with a 
standard reference, indium, and the reference used during the measurements 
was an empty copper DSC pan. 

Samples were pre-tested, by storing them in an oven at 80°C for two 
hours, and then immediately transferred to a _ dessicator. Pre-weighed 
samples. were. sealed. in. a...DSC. pan junder...pure Np. During the 
measurements, the samples were heated from ambient temperature up to 
600°C at the controlled heating rate of 40°C/min, and then weighed after 
completion of the reaction. 

A typical DSC result is shown in Fig. C.1 for an unretarded sample. 
A downward peak corresponds to an endothermic reaction, while an upward 
peak means an exothermic reaction. The total values of heat flow into (or 
out of) the samples during the reaction were automatically integrated by the 
DSC data recording and processing system. Fig. C.2 shows a comparison of 
the DSC results for pure and retarded samples. It is apparent that the 
unretarded sample went through endothermic reaction, even though there is 
a minor endothermic reaction at the beginning. 

The values of the heat of reaction per unit mass of volatiles were 


calculated by using the following formula: 


is 


r (C=1) 
1-Y 


where 


q = heat flow per unit mass of sample, 


obtained from the DSC integrated value. 


NG = mass fraction of char. 


The average values of the heat of reaction obtained from the DSC 


measurements are 0.23 x 10° J/kg for the unretarded samples (Yp = 0) and 


-0.10 x 10° J/kg for the retarded samples (Yp = 0.03). They are used in 


the numerical calculations of cellulose pyrolysis. 
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APPENDIX D 


CALIBRATION OF GAS CHROMATOGRAPH 


The purpose of the GC measurements 1s to obtain the mass fraction 
of noncombustible volatiles in the total volatiles for cellulose pyrolysis. The 
results of the GC measurements are given by the integrator as an area in 
counts for each component. Therefore a calibration factor is néeededwercr 
converting the area into the mass for each component. 

A bottle of calibration gas (Alltech Associates, Inc., Scott Specialty 
Gases, Can Mix 217) is used to obtain the calibration factors for CO and 
CO». (The factors for Oy and No» are also obtained at same time.) The 
bottle contains known volumetric percentages of CO, CO, and O> in an N>5 
base. The formula for calculating the mass of each component derived by 


using the perfect gas equation is 


P 
mj; = RT VeVi A (D-1) 
where m; = mass of species i 
P = pressure at which the calibration gas is sampled 
= | atm 
R = universal gas constant 


= 8314 J/kg mol.K 
T = temperature at which the calibration gas is sampled 


= ambient temperature 


eine 


= 
It 


: volume of the calibration gas sample 


nie 


M; molecular weight of species i 


x 


: volumetric fraction (i.e. mole fraction) of species i. 


The mass is divided by the corresponding "“counts" to give the 
calibration factors F; in units of ywg/count. Once the factor is known, the 
masses of CO, and CO for unknown samples can be calculated, simply by 
multiplying the total counts c; by the calibration factor F;. 

A mixture of laboratory air and pure N>y was used to calibrate for 
H»5O. During pyrolysis the H,0O contents for pure and fire retarded 
cellulose samples cover a wide range. Because the H5O content in ambient 
air is limited by the ambient temperature and relative humidity, it was not 
possible to get a sufficiently small HO content from the ambient air. 
Therefore, it was decided to mix laboratory air with pure N>. Thus the 
desired eect of H,O can be obtained by adjusting the air and N5 
volume ratio which includes the pure laboratory air case. The formula for 
calculating the mass of H,O in the mixture was obtained by using Van 
Wylen and Sonntag’s equation (1985), which gives the mass of Hj0O in air, 


and the ideal gas equation. It is as follows: 
where MH,O = mass of water vapor in the mixture of laboratory air and 


N) 


a ky. 


where 


ma = mass of laboratory air 


My = mass of the mixture of laboratory air and N>5 
Vim = volume of the mixture of laboratory air and N>5 
= 1 ml 
Rm = gas constant of the mixture of laboratory air and N5 


= 292 J/kgK (for 90% of N>) 


Tm = temperature of the mixture of laboratory air and N> 
Poat = Saturated vapor pressure of steam at T,, 
dm = relative humidity of the mixture of laboratory air and N> 


In the case of pure laboratory air, Ma=M,,, the Eqn. (D-2) becomes 


VA 
MH ,O = 0.622 ar oe Psat PA (D-3) 
A ‘A 
Va = volume of laboratory air 
= 1 ml | 
Ra = gas constant for air 
= 28 ats Kk SUS 
Ta = temperature of laboratory air 
Psat = saturated vapor pressure of steam at TA 


¢a = relative humidity of laboratory air 


According to the formulas, since Viet EA Ri and Ra are constants, 


and Peat is a function of T,,(T,) at one atmosphere, the mass of HO is a 


function of T,(T,) and $,,(¢,). Because po; is a very sensitive function 


cade 


of temperature, a thermometer with iore divisions was used to measure 
Tm(Ta). A relative humidity meter (Dikson Instant RH) was used to 
measure the $, (@a) values to 1 percent accuracy. chee meter was 
calibrated against 0% and 75% standard relative humidity capsules. 

After a number of calibration runs, the values of the mass of H50 
per unit volume were obtained as functions of a wide range of GC "counts". 
The calibration results of m vs. c of H5O are shown in Fig. 5.3. Since 
the calibration factor F; = mj;/c;, the results show that the calibration factor 
for H5O is not a constant. A linear fit was run through the calibration 
data. The equation for the straight line was used for data reduction. 

Once the calibration factors are known, the mass fractions of the 
noncombustible gaseous products of cellulose pyrolysis can be calculated as 


described below. 


The mass fraction of species i in total volatiles is defined as the 
ratio of mass flow rate of species i to the total volatiles production rate (as 


measured by the electronic balance during the experiment), Le. 


Virion (D-4) 


The basic assumption made for the calculation is that the mass 
fraction of species i in the total volatiles, as defined above, is a constant 
throughout the whole volatiles flow system (see Fig. 4.1 and Section 4.1) 


1.¢. 


Y= ¥ 


i i, Rotameter= Yi, sample loop (D-5) 


o15- 


There are the following relations. 


m; = Yi Rotameter ™Rotameter 


= Yisample loop ™Rotameter (D-6) 
Cj Fj 
Yisample loop ~ (D-7) 
Dy (c; ry) 
i 
MR otameter ~ VRotameter PRotameter (D-8) 


The value of Ppotameter iS taken as the value of PN, at 1 atm and 
the measured value of T because it is found that there is about 99% N> in 
the total volatiles during the cellulose pyrolysis. 


Therefore, the final formula is 


counts Bs 


VRotameter PN> (D-9) 


x (counts F;) 


The my data are obtained by using the same balance-computer 


ol 
system. The data acquisation and processing procedure is exactly the same 
aSsinathe hg measurements. 


At last, the data for m; and m,,) are combined to get the Y; 


values. 
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Table 2.1 


Reflectance of Sample Surface 


522- 


0.03 


0.04 


0.03 


0.03 


0.03 


TABLE 3.1 


Material Properties Used in Pyrolysis Calculations 


PS 


(kg/m?) 


550 


550 


Kg 0 
(w/m°C) 


0.1 


0.1 


E 


(J/kg mol) 


1.5x108 


1.5x108 


PC 
(kg/m?) 


55 


165 


KC @ 
(w/m°C) 


0.027 


0.080 


Cos 
(J/kg °C) 


1300 


1300 


h 


(W/m? °C) 


ees 


13.5 
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Coc 
(J/kg °C) 


900 


900 


* 
ke 


(w/m°C) 


4x1074 


4x1074 


0.96 


CoG 
(J/kg °C) 


1700 


1700 


h 


In 


(J/kg) 


-0.23x10° 


0.10x10° 


qab,s 
(kW/m2) 


40 


40 


TABLE 5.1 


Mass Fractions of Char from Cellulose Pyrloysis 


OS Yohar 
0 0.09 
0.005 0.19 
0.010 0.25 
0.015 0.27 
0.020 0.28 
0.025 0.29 
0.030 0.30 
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TABLE 5.2 


Elemental Composition of Ccllulose 


YC cell YH cell YO.cell 
0.44 0.062 0.49 
TABLE 5.3 


YR ~ YCchar 

0 0.80 
0.005 0.81 
0.010 0.75 
0.015 0.76 
0.020 Oa7 
0.025 0.77 
0.030 0.76 


Elemental Composition of Char 


YH char 


0.039 


0.041 


0.044 


0.042 


0.040 


0.041 


0.036 
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YO char 


0.16 


0.14 


0.18 


0.17 


0.15 


0.14 


0.14 


YNa,char 


Elemental Composition of Total Volatiles 


0.34 


0.32 


0.31 


0.30 


0.29 


TABLE 5.4 
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0.60 


0.61 


0.62 


0.64 


0.64 


TABLE 5.5 


Mass Fractions of Non-combustible Gases in Total Volatiles 


bas: mC. YHO > FT 

0 0.06 0.25 0.69 
0.005 0.16 0.37 0.47 
0.010 0.19 0.43 0.38 
0.015 0.21 ee 0.35 
0.020 0.22 0.44 0.34 
0.025 0.23 0.43 0.34 
0.030 0.24 0.41 0.35 


oT 


TABLE 5.6 


Elemental Composition of Combustible Volatiles and Stoichiometric Ratio 


YR Tel CV x CV Yo, CV Vg 

0 0.57 0.050 0.38 1.6 
0.005 0.66 0.050 0.29 is 
0.010 0.75 0.050 0.20 2.2 
0.015 0.75 0.060 0.19 2A) 
0.020 0.73 0.060 Ue ie 
0.025 0.68 0.060 0.26 51 
0.030 0.63 0.070 0129 2.0 
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TABLE 5.7 


Major Products of Ccllulose Pyrolysis 


0.28 


0.30 
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0.47 


0.46 


Y 


com,vol 


0.24 


TABLE 6.1 


Material Properties Used in Diffusion Flame Calculations 


0.61 


0.46 


0.40 


0.38 


0.38 


O77 


0.34 


0.35 
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Schematic drawing of the apparatus used 
in heat of gasification measurements. 
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Figure 2.2 Schematic drawing (top view) of set-up 
for cellulose pyrolysis measurements. 
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Figure 2.4 Vertical variation of absorbed heat flux 


over the sample height. 
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Figure 2.5 Voltage outputs of thermocouple and 
pyrometer vs. time. 
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Figure 2.6 Calibration of surface temperature vs. pyrometer 
Output voltage, the data is fitted by the polynomial 
T = -116 + 641V - 157V2 + 14.6V2, 
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Figure 2.7(a) An example of computer output of 
surface temperature vs. time; Yp=0; 


Gab,s=40 kW/m? 
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Figure 2.7(b) An example of computer output 
mass loss rate vs. time; Yp=0; 


da bs =40KW / nr. 
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Figure 2.7(c) An example of computer output of 
heat of gasification vs. time; Yp=0; 


ap, 5=40 kW/m?. 
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Figure 2.8 Surface temperature vs. time for different 
absorbed heat fluxcs; Yp=0. 
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Figure 2.9 Mass loss rate vs. time for different 
absorbed heat fluxes: YRp=0. 
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Figure 2.10 Heat of gasification vs. time for different 
absorbed heat fluxes; Yp=0. 
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Figure 2.11 Surface temperature vs. time for pure and 


retarded cellulose samples; q'y .=40 kW/m2, 
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Figure 2.12 Mass loss rate vs. time for pure and retarded 


a ; p) 
cellulose samples; q4b,s749 kW/m°~. 


-105- 


3.0 


2.5 
2.0 
oO 
aw 
= 
=e) 
o 
(S) 
= 15 
x 
oOo 
of exes 
Y 
1.0 R 
9) 
0.5 
0.03 
O 60 120 180 240 300 


t(s) 


Figure 2.13 Heat of gasification vs. time for 
pure and retarded ccllulose samples; 


on" 


D 
Vab,s=49 kW/ntr. 
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Figure 2.14 Surface temperature vs. retardant concentration; 


Gab.s=40 kW/m; t=200s. 
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Figure 2.15 Mass loss rate vs. retardant concentration; 


q" =40 kW/m2: t=200s. 
ab,s 
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Figure 2.16 Heat of gasification vs. retardant concentration; 


q’-=40 kW/m?; t=200s. 
ab,s 
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Figure: 3.1 


Control volume used in pyrolysis analysis 


Olea charring fuel: 
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Figure 3.2 Theoretical surface temperature vs. time 


for pure and retarded cellulose. 
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Figure 3.3. Theorctical mass loss rate vs. time 


for pure and retarded cellulose. 
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Figure 3.4 Theoretical heat of gasification vs. 


for pure and retarded ccllulose. 
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Figure 3.5 Experimental and theoretical mass loss rate 


vs. time for pure ccllulose; Yp=0. 
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Figure 3.6 Experimental and theoretical mass loss rate 


vs. time for retarded cellulose; Yp=0.03. 
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Figure 4.1 Schematic drawing of apparatus used in 
measurements of gaseous products. 
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Figure 4.3 An example of output 
Y p=0.03; t=180s. 


-118- 


B35 32 
ear 3 
yt soo 
gos Soz 
emia 

b> a= =, 
er si, 
ei Sree Sis 
ye 75 
i <- = 
aos ‘ 
SR ere 
a 

ws 

Bere od 
Seige ih 
ey Se 
g25 32 
Se) or 
le sy eee 


of GC analysis; 


Ol6l-ZEE “ON LWVd HSWIS-NIMdSc 


0 90 180 270 


Figure 4.4 Mass fraction of CO5, H,0 and CO 


vs, time for pure and retarded 


cellulose samples; Gab.s=40 kW/m2 
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Figure 4.5 Mass fraction of CO, HyO and CO 


vs. retarded concentration; 


Gy p,s=40 kW/m* ; t=180s. 
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Figure 5.1 Effects of retardant concentration on major 
products of cellulose pyrolysis. 
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Figure A.2 Schematic drawing of the apparatus 


for making cellulose samples. 
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Figure C.2 Comparison of DSC results for pure and 


retarded cellulose samples. 
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Figure D.1 Calibration curve of water mass vs. integrator 


“count. 
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